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ABSTRACT 


Ten  algorithms  pertaining  to  underwater  acoustics, 
target  motion  analysis,  P-3  mission  planning,  flight  crew 
management,  and  naval  gunfire  support  conversions  are  presented 
along  with  programs  for  Hewlett-Packard  HP-67  and  HP-97  pro- 
grammable calculators. 


I . INTRODUCTION 


k 


This  report  contains  a collection  of  programs  which  were 
submitted  by  officers  in  partial  fulfillment  of  the  requirements 
of  the  course  Tactical  Design  and  Analysis  (OA  4658)  conductec 
at  the  Naval  Postgraduate  School  during  the  period  of  October 
through  December  1978. 

All  programs  were  listed  using  an  HP-97  with  HP-97  key 
codes.  The  corresponding  HP-67  key  codes  may  be  found  on  pages 
324  through  331  of  the  "HP-67  Owner's  Handbook  and  Programming 
Guide. " 


II. 


ACTIVE  SONAR  ACQUISITION  by  Mr.  R.  F.  Fish  and  LT  M.  H.  Trent 


A.  Problem  Statement 

A sonar  at  a depth  (SD)  has  the  possiblity  of  detecting 
a target  at  a depth  TD  at  a slant  range  r' . Detection  can  only 
occur  if  the  target  lies  within  the  beam  pattern  and  the  signal 
excess  is  at  least  equal  to  the  detection  threshold.  Whether 
or  not  the  system  is  noise  or  reverberation  limited  depends  on 
the  geometry  and  doppler  frequency  shift.  The  problem  is  to 
determine  the  acquisition  range  of  the  sonar  with  various  geo- 
metries and  acoustic  parameters. 


B.  Operational  Analysis 

The  analysis  uses  a 0 dB  detection  threshold  because 
of  the  limited  number  of  storage  registers  (26)  and  program 
steps  (224)  in  the  calculator. 

In  using  the  program  it  should  be  noted  that  calculations 
do  not  include  the  effect  of  shadow  zones.  Acquisition  ranges 
computed  must  be  considered  with  this  in  mind.  In  addition,  once 
the  target  and  surface  signals  are  outside  the  beam  pattern 
(+  3db)  they  are  assumed  to  abruptly  disappear  respectively. 

The  analysis  also  assumes  the  water  is  deep  with  no  bottom  effects. 

Considering  the  above  caveats  the  source  and  target  can 
be  placed  as  desired  in  the  medium  and  the  appropriate  sonar 
equation  parameters  will  be  computed.  In  doing  so,  several  tests 
will  be  made  to  determine  which  equations  will  be  used. 


PRECEDING  FiflE 


The  program  will  terminate  before  acquisition  if  one 
of  the  following  occurs: 


Slant  range 
r'  < 0, 

Angle-to-target 

0^  > <t>/2  + y,  or  ©i  > 4 >/2  - y (if  target  is  below  source)  , 

Angle- to-surface  at  r' 

02  > 4>/2  + y, 


where 

y = pitch  angle  range . 


These  terminations  and  signal  excess  SE  _>  0 will  finish 
with  a "1"  printed  as  an  output  at  the  end  of  the  calculations. 

After  calculation  of  the  surface  reverberation  level, 

RL  , there  is  a program  stop  where  the  appropriate  correction 

S 

can  be  input  for  off-axis  transmission  and  reception.  The  same 
event  occurs  when  SF  is  calculated  so  that  the  sonar  equation 
can  be  corrected. 

The  output  listing  is  as  follows: 

Doppler  No  Doppler 


r’ 

91 

e2 

TL 

NL 

SE 


s 


r' 

61 

02 

TL 

RL  displays  but  not  printed 

S 

Total  of  RL  + RL  + NL  combined  and  SE 
S V s 


4 


c. 


Computational  Algorithm 


1 . Input 

— Sound  speed,  c (m/sec) 

— Listening  time  between  transmit  pulses,  t (sec) 

--Source  depth,  SD  (m) 

— Target  depth,  TD  (m) 

--Horizontal  (and  vertical)  half-beam  width,  4>/2  (degrees) 
— Sonar  pitch  angle,  y (degrees) 

--Mixed  layer  depth,  D (meters) 

— Absorption  coefficient,  a (dB/meter) 

— Frequency,  f (kHz) 

— Sea  state,  S.S. 

— Column  scattering  coefficient,  Sc 
— Constant,  10 

--Sonar  self-noise  level,  NLg  (dB) 

--Target  strength,  TS  (dB) 

— Range  decrement  (meters) 

— (Pulse  length,  x (sec)  * 0 (radians)]  2 
--Surface  scattering  coefficient,  Sg  (dB) 

— Sonar  source  level,  SL  (dB) 

2.  For  doppler  set  Flag  0;  for  no  doppler  clear  Flag  0. 

3.  Output 

r',  0^,  0J#  TL,  RLg,  RLv>  total  level,  SE.  For  doppler, 

RL  , RL  , and  total  level  are  included  in  NL  . 
s v s 
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D.  HP-67/97  Calculator  Program 


1.  User  Instructions 


Instruction 


Enter  program  card 

Enter  data  in  primary  stage: 
Sound  speed  (m/sec) 

Listening  time  (sec) 

Source  depth  (m) 

Target  depth  (m) 

Half-beam  width  (deg) 

Sonar  pitch  angle  (deg) 

Mixed  layer  depth  (m) 

Absorption  coefficient  (dB/m) 
Frequency  (kHz) 

Sea  state 
Source  level  (db) 

Enter  data  in  secondary  storage:' 
Column  scatterinq  coefficient 
Constant 

Sonar  self-noise  level 
Target  strength 
Range  decrement  (m) 

(t  (sec)  x 4>(radians)  l/2 

Surface  scattering  coefficient 

Primary/Secondary  exchange* 

Doppler  or 
No  doppler 

Start  computations 

Printed  output: 

Slant  range 
Angle -to- target 
Angle- to-sur face  at  range  r' 
Transmission  Loss 
If  Flag  2 is  set,  self-noise 
is  printed. 


| Key(s) 

Output 

STO 

0 

STO 

1 

STO 

3 

STO 

4 

STO 

7 

STO 

8 

STO 

A 

STO 

B 

STO 

C 

STO 

D 

h STI 

P ? 

S 

STO 

0 

STO 

2 

STO 

3 

STO 

4 

STO 

5 1 

STO 

6 i 

STO 

8 

P Z 

S 

SF  0 

doppler 

CF  0 

no  doppler 

A 

See  Step  7 

r' 

*1 

02 

TL 

NL 

s 

Step 

Instruction 

Input 

Key (s) 

Output 

8. 

If  Flag  2 is  set,  go  to  Step  9. 
Otherwise : 

a. 

Display  surface  reverberation 

RL 

level 

s 

b. 

Enter  two-way  beam  pattern 

correction  in  db  (0  if  no 
correction) 

Correction 

R/S 

Total 

Print  total  corrected  level. 

RL 

s 

9. 

Stop  and  display  SE.  To 

SE 

change  range  decrement, 
execute  Step  a.  Other- 
wise go  to  Step  b. 

a. 

Key  in  new  range  decrement 

r.  d. 

STO  5 

R 

SE 

b. 

Enter  two-way  beam  pattern 

correction  in  db  (0  if 
no  correction) 

Correction 

R/S 

Corrected 

Display  corrected  SE 

SE 

10a. 

If  SE  < 0,  execution  con- 

tinues  from  Step  7. 

b. 

If  SE  0,  termination  occurs 

1.00 

The  primary  and  secondary  registers  must  be  exchanged  before  (Step  3) 
and  after  (Step  4)  entering  data  into  the  secondary  storage  registers 


Sample  Problems 


i.  Sonar  and  tarqet  are  in  the  mixed  layer  in  a "doppler" 
situation,  go  that  acquisition  ia  noise  limited  by  NL 

s 

(Figure  1) . 


Source  depth  = JO  meters 
Tarqet  depth  = (>0  meters 
Mixed  Layer  depth  = 75  meters 


FIGURE  1 . 


Geometry  of  Sample  Problem  1. 


Input . 


R0  : 

1500  m/sec 

SO  : 

-50  dB 

R1 : 

6 . 7 sec 

S2  : 

10 

R3s 

30m 

S 3 : 

65  dB  re  In Pa 

R4  : 

60m 

S4  : 

10  dB 

R7 : 

10” 

S5  : 

500  m (initial 

R8 : 

5 * 

S6  : 

. 1 7 "see 

RA: 

75m 

S8 : 

-30  dB 

RB: 

.00328  dB/m 

RC: 

25  kHz 

RD: 

2 

RI  : 

22  7 dB  re  1 yiPa  0 1 m 

Set  Flag  0.  No  corrections  are  made  to  SF . 
Output 


ramie  decrement  ) 


The  results  are  shown  on  the  Sample  Problem  output . 
At  acquisition 

r'  * 3650  meters 
6j  « .47 
02  = .47 

TL  = 85.1  dn  (layer) 

NL^  = 65  dB  re  lpPa 
SF.  =>  .21  dB 


1 


SF0 

Set  for  doppler 

G£Sn 

Start 

5025.0*3 

**»■ 

r' 

0.34 

**♦ 

el 

0.34 

*** 

02 

53.  76 

444 

TL 

05.00 

444 

RLS 

0.00 

R.£ 

Correction  to  SE 

-15.52 

4M 

Corrected  SE 

4525.00 

0.30 

*** 

0.38 

444 

50.55 

»*» 

05. 00 

*4* 

0.00 

R 2 

-5.50 

***• 

4025.00 

**» 

0.43 

*** 

0.43 

8S.05 

*** 

85.00 

jM* 

300.00 

C^e-  "I 

Change  range  decrement 

Ri  J 

to  300  meters 

0.00 

R/£ 

-4.17 

44* 

3725. 08 

44* 

0.46' 

444 

8.4b' 

*** 

86.34 

444 

65. 00 

44* 

50.09 

£705  ] 

Change  range  decrement 

R*  J 

0.00 

R.  £ 

-0.  bl 

*44 

3875.00 

4*4 

0.47 

*** 

0.47 

8b'.  04 

**4 

85.00 

44  4 

25.00 

£705  1 
R1  J 

Change  range  decrement 

0.00 

R-  £ 

-0.05 

444 

3650. 00 

$4-4 

r’ 

0.47 

44  4 

01 

0.47 

444 

02 

85.56 

44  4 

TL  (layer) 

65.00 

444 

NLS 

0. 00 

R.  £ 

0.21 

**> 

SE  = 0.21 

1.00 


SE  = 0.21 
Terminate 


Sample  Problems  2 and  3.  The  sonar  and  target  configuration 
are  shown  in  Figure  2. 


Surface 
Mixed  Layer 


Source  depth  = 600  meters 
Target  depth  = 600  meters 
Mixed  Layer  depth  = 75  meters 


FIGURE  2.  Geometry  of  Sample  Problems  2 and  3. 


Input  for  Problems  2 and  3: 


Same  as 

for  Problem  1 except 

R3 : 

600 

meters 

Flag  0: 

set  for  problem  2; 

R4  : 

600 

meters 

clear  for  problem  3 

R8: 

0° 

0 entered  as  correction  factors  to  RL  and  SE. 

s 

Output  for  Problem  2: 
r'  = 4125  meters 


92  = 8.36  ° 

SE  = .32  dB 

Sonar  acquired  the  target  at  abou«_  4125  meters. 

Output  for  Problem  3: 

r'  = 3025  meters 

ex  = o° 

e2  = ii.44° 

SE  = 0.85  dB 

Sonar  acquired  at  3025  meters. 


13 


SAMPLE  PROBLEM  2.  Output 


14 


CF7 

No  doppler 

t»6£  H 

Start 

5025. Oo 

*»  * 

r' 

0.00 

Ml 

0.  60 

**» 

02 

.op.  50 

Hi 

TL 

77.67 

0 » ♦ 

RLS 

0.00 

£ £■ 

No  correction 

» i • w'Y 

• Of 

Total  RI 

0.00 

K £ 

s 

No  correction 

-2j.se 

♦ 99 

Corrected  SE 

4525.00 

*M 

0.00 

»** 

7. 62 

Ml 

27.35 

HI 

St.  71 

0.00 

R.6 

81.05 

m 

0.00 

* S 

-20.76 

IM 

4025.00 

HI 

0.00 

*»* 

8.57 

Ml 

85.  S0 

m 

86.  52 

• 0 0 

0.00 

R 2 

88.53 

m 

0.00 

R'S 

-20.13 

Ml 

7525.00 

0.08 

*»i 

3.80 

in 

62.51 

Ml 

31.53 

Ml 

0.00 

R S 

31.58 

999 

0.00 

-18.53 

in 

3025. 00 

99  0 

r' 

0.00 

in 

01 

11.44 

990 

02 

73. 54 

mi 

77. 06 

»»? 

0. 65 

in 

0.00 

s 2 

0. 85 

*•» 

SE 

1.00 

Termination 

SAMPLE  PROBLEM  3.  Output 


3.  Program  Storage  Allocations  and  Program  Listings 


Registers 

R0: 

c (m/sec) 

SO: 

Sc  (dB) 

Rl : 

t (sec)  and 

SI : 

RLv  (dB) 

SL  - 2TL  + 10  log  r' 

S2 : 

10 

R2: 

R and  r' 

max 

S3: 

NL  (dB  . _ ) 
relyPa 

R3 : 

SD  (m) 

S4  : 

TS  (dB) 

R4  : 

TD  (m) 

S5: 

Range  decrement  (m) 

R5: 

6^  (deg) 

(500  m default) 

R6  : 

02  (deg) 

S6  : 

t (sec)  x <p/2  (AD) 

R7: 

<f»/2  (deg) 

S7: 

10  log(4|tc/2) 

R8 : 

Y (deg) 

S8  : 

Ss  (dB) 

R9 : 

<*/2  + y (deg) 

S9  : 

RL  (dB) 
s 

RA: 

D (meters) 

RB: 

a (dB/meter) 

RC: 

f (kHz) 

RD: 

S.S. 

RE: 

TL  (dB) 

RI : 

SL  (dBreluPa  * lm) 

or  consistent  with  NL 


Initial  Flag  Status  and  Use: 

0:  ON  for  doppler,  1,  2,  3:  OFF,  unused 

OFF  for  no  doppler 

16 
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E.  Computational  Analysis 

The  active  sonar  equation  is 


where 

SL 

TL 

(NL-DI ) 

RL 

DT 

TS 

and  NL 


SL  - 2TL  + TS  - (NL  - DI ) _>  DT  , 

NL 


RL 


source  level  for  the  sonar  ^BrelviPa  ® ' 

transmission  loss  (dB) , 

ambient  noise  term  which  is  neglected  due  to 
a higher  NLg  or  RL, 
reverberation  level  (dB) 

system  detection  threshold  ( OdB  assumed), 

target  strength  (dB) , 

Self-noise  (dB  , _,  ) . 

relpPa 


The  only  terms  not  known  in  the  equation,  TL  and  RLg, 
are  calculated  at  various  ranges  (decrements)  until  the  signal 
excess  (SE)  >_  DT(0)  . 

The  TL  is  calculated  for  two  conditions: 
a.  When  both  source  and  target  are  in  the  layer  (Reference  1) 


TL  = 10  log  r.  + 10  log  r'  + ar1  + , 

^ t ^ rs 


where 
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105 


is  the  transition  range  (meters) , 


r 


t 


a = absorption  (dB/meters) , 

b = 1.04  x SS  x /f  bounce  factor  (dB /bounce)  valid 

between  (3-25  kHz) (3-14  dB/bounce) , 

rs  = 840  /D  , 

zg  = larger  of  source  or  receiver  depths  (meters)  , 
and  D = layer  thickness  (meters). 


b.  When  both  source  and  target  are  not  in  the  layer  the  TL  is 

TL  = 20  log  r'  + ar'  (for  r'  < rfc  also  when  in  layer). 

After  using  the  proper  TL  formula  it  must  be  decided 
whether  or  not  there  is  sufficient  doppler  to  be  able  to  dis- 
regard reverberation  (i.e.,  RLg  and  RLv) . If  there  is  enough 
doppler  then  the  NLg  term  dominates  and  the  sonar  equation  can 
be  solved.  If  there  is  no  doppler  the  appropriate  reverberation 
must  be  considered  and  combined  with  NLg.  Then  the  sonar  equation 
can  be  solved.  By  successive  decrements  of  r' , there  may  be  a 
point  where  SE  ^ DT  and  thus  detection  has  occurred. 

The  reverberation  equations  are  (Reference  2) , 

RLg  = SL  - 2TL  + 10  log  r'  + Sg  + 10  log  (-^-)  , 

where 
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...ii 


S,  = surface  scattering  parameter  (dB)  for  the  particular 
conditions  (wind  speed,  grazing  angle) , 
ij)  = sonar  horizontal  beam  width  (radians)  , 
c = wave  propagation  time 
and  x ■ transmit  pulse  width  (seconds) ; 

and 

RLv  = SL  - 2TL  + 10  log  r'  + + 10  log  (^y-M  , 

whe  re 

Sc  = column  scattering  coefficient  (dB)  for  the  particular 
environmental  conditions. 


At  each  range  decrement  and  0^  are  calculated  to 

determine  if  they  are  inside  the  beam  pattern.  If  is  not, 

acquisition  cannot  occur.  If  0,  is  not,  then  RL  is  not 
important.  The  formulae  for  these  quantities  are  (Figure  3): 


Pj  = sxn 


-1 


SD  - TD  , 

1 and 


sin 


-1  SD 


where 

SD  = source  depth, 
and  TD  = target  depth. 


In  addition,  depending  on  the  values  of  these  angles, 

corrections  can  be  made  to  the  RL  and  sonar  equation  to  com- 

s 

pensate  for  off-axis  (beam  pattern)  transmission  and  reception. 
To  get  an  initial  value  of  r'  = R the  equation 

Hid  X 

R^ax  = ct/2  is  used  where  t = the  sonar  "listening"  time  or 
time  between  successive  pulse  transmission. 


2 2 
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III.  THREE  ENGINE  AVAILABLE  RANGE  REMAINING:  P-3(B)  AIRCRAFT 


CONFIGURATION  11 B"  by  LT  R.  J.  Knight 

A.  Problem  Statement 

Aircraft  total  fuel  remaining,  outside  air  temperature, 
and  aircraft  altitude  data  are  available.  Determine  the  air- 
craft's available  range  remaining. 

This  program  allows  a pilot  or  copilot  to  rapidly  and 
efficiently  provide  a quick  estimate  of  available  range  remaining 
in  an  emergency  situation  (three  engine  flight) . 


B.  Operational  Analysis 

The  aircraft's  available  range  remaining  can  be  extracted 
from  the  table  listed  on  pages  12-189  of  the  P-3(B)  aircraft 
NATOPS  manual. 


C .  Computational  Algorithm 

1.  Input  fuel  remaining  (pounds). 

2.  Input  outside  temperature  (°C). 

3.  Input  altitude  (feet) 

4.  Calculate  the  three-engine  available  range  remaining. 
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D.  HP-67/97  Calculator  Program 
1.  User  Instructions 


Step 

Instruction 

Input 

Key (s) 

Output 

1. 

Load  program  magnetic  card 
side  #1  and  side  #2 

0.000000000 

2. 

Input  fuel  remaining,* 
press  enter 

pounds 

ENT 

pounds 

3. 

Input  outside  air  temperature 
for  specific  altitude 
press  enter 

°C 

ENT 

°C 

4. 

Input  altitude 

feet 

feet 

5. 

Press  A to  calculate  the  three 
engine  available  range  remaining 

A 

range  in  NM 

Note:  Fuel  remaining  must  be  entered  as  #10,000,  #20,000,  #30,000,  #40,000 
or  #50,000.  "Error"  will  display  otherwise. 


2 . Sample  Problem 

Calculate  the  three-engine  available  range  remaining 
if  the  remaining  fuel  is  10,000  pounds,  the  outside  temperature 
is  11 °C,  and  the  altitude  is  2000  feet.  (ANSWER:  295  n.mi.) 


1660C.  ENT! 
11.  ENT! 

2600.  GS£n 
295.  **< 
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3 . Program  Storage  Allocation,  Permanent  Data,  and 
Program  Listing 
Registers : 


RO: 

a0 

SO: 

A: 

Fuel 

Rl: 

b0 

SI : 

B: 

Altitude 

R2: 

al 

S2 : 

C: 

Temperature 

R3: 

bl 

S3: 

D: 

Temperature  deviation 

R4: 

a2 

S4  : 

E: 

Uncorrected  range. 

R5  : 

b2 

S5  : 

R6 : 

a3 

S6: 

R7: 

b3 

S7: 

R8: 

a4 

S8 : 

R9  : 

b „ 

S9  : 

Flags : OFF,  unused. 


User  Control  Keys: 
A:  Compute 

B: 

C: 

D: 

E: 


a . 


b 

c 

d 

e 
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Permanent  Data 


The  following  permanent  data  are  stored  in  the  primary 
storage  registers  RO  through  R9 . 


-343l\2..?t>$lt*4  t 
122.  99t?3?6 
-35346.62(33 
43.1 8€H?69  3 

-3PS53.  e«4;093  4 

2?.m4H?i  r 
-3f^t'8.i5(?5tl  c* 
2*. 1S31S143  7 

-33432.66271  6 

It. 59293261  $ 


i 


i 


I 


Program  Listin 


Cl 

-45 

16-43 

— 4 
O 1 

2*  0: 

1 

61 

- 

-45 

,V«0? 

16-42 

ST  04 

22  64 

1 

61 

- 

-45 

,V*0  ■> 

16-45 

£705 

\'r> 

0 

+ 

ilBLI 
RC  LB 
KCL0 


-41 

00 

-24 
21  t. 

3c  12 
36  00 
-45 


656 

i 

-24 

03? 

STOE 

35  15 

646 

STO* 

22 

16  15 

641 

*LBL  2 

21  62 

645 

RCLB 

36  12 

647 

Rl\2 

36  65 

044 

- 

-45 

045 

RSL3 

36  63 

046 

T 

-54 

047 

57£lf 

35  15 

046 

C76* 

22 

16  15 

! 04? 

t.e  L3 

21  63 

050 

RClB 

36  12 

651 

SCI  4 

36  64 

655 

- 

-45 

655 

p:l5 

36  65 

654 

+ 

-i4 

i 655 

ST£I£ 

35  15 

056 

fT0» 

22 

16  15 

657 

♦ LF14 

51  04 

START 

STO  ALTITUDE 

SIX')  TEMT 

STO  FUEL 

TEST  FUEL 
- 10,000  lbs? 


TEST  FUEL 
- 20,000  lbs? 


TEST  FUEL 
= 30,000  lbs? 


TEST  FUEL 
*=  40,000  lbs 


TEST  FUEL 
- 50,000  lbs 


TEST  FUEL 
> 50,000  lbs 


COMPUTE  RANOE 
W/O  TEMP  CORRECTION 
F «=  10,000  lbs 


COMPUTE  RANCIE 
W/O  TEMP  CORRECTION 
F - 20,000  lbs 


COMPUTE  RAN  OF. 

W/O  TEMP  CORRECTION 
F - 30,000  lbs 


. Mathematical  Analysis 


A linear  curve  fit  was  performed  using  the  HP-67/97 
standard  pack  SD-03A  program.  Five  "fits"  were  preformed. 
For  a constant  fuel  weight,  X represented  the  range  and  Y 
represented  the  altitude.  Resulting  outputs  provided  the 
following: 


Fuel  (pounds) 


10,000 

20,000 

30.000 

40.000 

50.000 


,998615613 


998739204 


999756772 


999823326 


999802631 


-34,302.26904 

-35,846,62839 

-36,857.62093 

-38,068.60561 

-39,482.60271 


122.9065370 


43.18615709 

27.00641670 


20.18918143 


16. 50293201 


For  a constant  fuel  the  range  X could  be  obtained  as  follows 


Y - a 


„ altitude  - intercept  (constant  fuel) 

Range  slope  (con  si  ant'  Fuel) 


Temperature  correction: 

increase  range  1%  per  5°C  above  standard 
decrease  range  1%  per  5°C  belcw  standard 
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Based  on  a lapse  rate  of  -2°/1000ft  the  standard  temperature 
for  a specified  altitude  is  obtained  by: 


-2 


altitude 

IMS 


+ 15  = standard  temperature  (for  altitude) 


Subtracting  the  input  temperature  from  the  computed  standard 
temperature  yielded  the  difference  in  °C  from  standard. 

A correction  factor  could  be  obtained  for  each  difference 
times  .002/degree  times  uncorrected  range.  The  range  could  then 
be  computed  by  summing  the  correction  factor  and  the  uncorrected 
range  value. 


I 
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F.  Program  Flowchart 


YES 


! CORRECTION-O 
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IV.  MISSION  PLANNING;  P-3(B)  AIRCRAFT  CONFIGURATION  "B" 

LT  M.  D.  Thomas. 

A.  Problem  Statement 

In  order  to  carry  out  the  various  operational  missions 
assigned  the  P3-B  Aircraft,  effective  utilization  of  the  plat- 
form is  essential.  All  aspects  of  the  mission  must  be  carefully 
planned.  Fuel  planning  directly  influences  endurance  and  the 
effectiveness  of  the  mission.  The  NATOPS  manual  provides  charts 
for  this  purpose.  Two  vital  charts  for  planning  are: 

1.  four  engine  maximum  range  operating  table?  used  in  pro- 
ceeding to  the  operational  area. 

2.  three  engine  loiter  operating  table;  used  while  onstation 
for  minimum  fuel  consumption. 

The  pilot  or  flight  engineer  enters  with  the  aircraft's  altitude 
and  gross  weight  and  finds  the  correct  indicated  airspeed  (IAS) 
to  fly. 

This  program  is  a user's  program  in  that  it  translates 
these  two  charts  onto  an  HP-67/97  magnetic  card  and  allows 
calculation  of  IAS  without  the  charts.  Most  missions  are  flown 
in  configuration  'B!  therefore  the  program  presented  here  is 
for  that  case. 

B.  Operational  Analysis 


None . 


C.  Computational  Algorithm 


1. 


! 

[ 


2. 


Enter  altitude  and  gross  weight  in  packed  form:  AAAAA. WWWW 
where  AAAAA  denotes  the  altitude  in  feet,  and  WWWW  is  the 
gross  weight  divided  by  100,000.  The  leading  zeroes,  if 
any,  in  the  value  of  WWW  must  be  entered.  For  example, 

18,000  feet  and  76,500  pounds  are  entered  as  AAAAA  = 18,000 
and  WWWW  = 0765,  that  is  18,000.0765. 

Compute  the  maximum  range  IAS  or  the  three-engine  loiter  IAS. 


D.  Hf-67/97  Calculator  Program 
1 . User  Instructions 


step 

Instruction 

Input 

Key (s) 

Output 

1. 

Enter  program  card 

2. 

Enter  altitude  and  gross  weight. 

G.W.  must  be  at  least  a three 
digit  number.  Compute  four  engine 
max  range  IAS . 

ALT. GW 

A 

IAS 

(max  range) 

3a. 

Enter  altitude  and  gross  weight. 

G.W.  must  be  at  least  a three 
digit  number.  Compute  three 
engine  loiter  range  IAS 

ALT. GW 

B 

IAS 

(Loiter) 

♦ 

1 


b. 


Optional : Compute  the  four  engine 

max  range  IAS  without  re-entering 
the  altitude  and  weight 


R/S 


IAS 

(max  temp) 


B 

Loiter  IAS 


2 .  Sample  Problems 

1.  130,000  lbs  at  18,000  ft 

18000.130  


2.  86,000  lbs  at  6,000  ft 

6000.086*  241 

3.  76,500  lbs  at  10,000  ft 

10000.0765  232 

4.  50,000  lbs  at  3,000  ft 

3000.050  Error 

5.  130,000  lbs  at  30,000  ft 

30000.130  Error 


A 

Max  Range  IAS 

252 


Gross  weight  must  be  at  least  a three  digit  number  (IE) 

130,000  .130 

76,000  .076 

82,500  .0825 


3.  Program  Storage  Allocation  and  Listing 


Registers 

R0:  altitude 

SO : 

RA: 

Rl: 

max  range  constant 

SI: 

R2 : 

loiter  airspeeds 

S2: 

• 

• 

• 

R9 : 

S9: 

RE: 

RI: 

35 


220 

175 

165 

Error 

Error 
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Initial  Flag  Status  and  Use: 
0:  OFF,  unused 

1.  OFF,  unused 


2:  OFF,  unused 
3:  OFF,  unused. 


User  Control  Keys: 


A: 

Compute 

four  engine  IAS 

a : 

B: 

Compute 

three  engine  IAS 

b : 

C: 

c s 

D: 

d: 

E: 

e : 

36 


- - ^ 


> V) 

tn 

&,  H 

o 

8 a 

a H 

H 

H tO 

to 

< s 

w 

u o 

p 

u 

1-1 

X 

< 

6 to 

g 

£ s 

M 

S w 

H 

S 

M 

to  w 

Q 

2 

to  S 

9 ^ 

8 

U 

2 2 

O M 

bu 

a.  in 

o 

8 § 

H 

to 

< in 

° § 

h 8 

J 

< 

0 to 

z 

m t 

o 

w g 

M 

s w 

t- 

s 

M 

to  w 

Q 

to  S 

2 

o o 

8 

a.  z 

u 

tj  M 

b » wo  *io  » « ro  io  ii  > w~»  ro  f i sr  *>4  *r  **■  wo 


•a  io  *a  *aa  ‘a 

i 

•O  IO  VD  S© 

-4  »\j  ro  m 


*'•>  *0  't  »VJ  •»  '*4  4*4  * 4 ‘Vi  IO  • » >45  'i*  '-O  IO  *a 

I I I • I II  I 

‘O  t\J  wo 

<M  ^ c\4  *«r 


' 4 Ii  I • 4 Ii  > lo  -f  Ii  1 II  I 


in  •"  *o  *x>  — • :<  wo 

‘U  «4  -J  U4 
41  40  <«>  U U4 
> . • Xi  lV  c* 


U“J  -E  ^ C « 4 • CO  IO  — • -♦  WO  «>4 

X I-  -a  >-  O I U-  I 

O IV  ui  Nl  I-  I—  K- 

— J >C  • -X*  •/»  Oi 


•X»  io  •© 
I 

*■•  '©  »vj 

i\i  *-4  C*4 


wo  x wo  c — . 

u I * — J v-  I*J 

►—  U:  io  *41  I— 
•l»  --I  ci» 


>•4  *'a  “4  « wo  r,-j 

• o ‘-ij  *:a  *3a  ro  •■©  •’a 

I I 

in  m 


—»  ♦ m i\j 
i iv  i 

►-  i— 

wi  u> 


05  0%  «Xi  i>4  ro  in  *-ai  i - o,  <5  •>«  *>j  ro  wo  '£•  r to  «J.  *x*  *-«  cvj  ro  o y , • ^ f ,y,  i»,  «i.,  rVJ  r*. 

O’  ^ ^ 7 •»-  m -r  -t  -t  *1-1  *1-1  u J u-i  in  Ii-.  *1- . u . m ir>  'ii  u.  'u  ‘a  'o  <■*>  •*:.  •*-,  u.  ‘L,  r r.  r * r . rC 
« ® •*»  «x»  ® *x*  »x»  •x*  cx*  «x>  ct*  taa  tii  *Xi  «x>  «4j  vi»  *x»  Ch.«  *3a  »i»  *i»  «3a  »aa  csa  «xi  'X»  *3>  «aa  ct»  cl»  *xi  *x»  ti<  *x  ct»  co 


§ s 

K 

X H 

U in 

a a 

u 


>■4  VJt 

• -*  ».  a 

* i 

*N4 

-r 

IO 

• a 
•o 

* 4 -r*  t o 

IO  ^ VnJ 

1*4 

• a 

l.  } 
IO 

•M 

•X) 

• 4 

•o 

1 - 
•*a 

» -4 

•'a 

4 4 

v v4 
•O 

• '4 
•S 

•a 
• a 

• 1 ' 4 

•o  io 

- 4 

•0) 

' i 

••a 

r - 
•3d 

• i 

•a 

U5 

•a 

'i 

•*a 

• V4 

* ^ 

»-4 

• o 
•'a 

'i 

vs4 

- 4 r . 
•‘a  O 

» 4 

ro 

» 1 
•a 

»-4 

•u 

V» 

*S4 

•o 

*•*  to 

'O 

in 

SO 

wo 

in 

••i 

•■■a 

IO 

*o 

— « r»> 

»\j  *o 

1 

•-4 

r«* 

»-4 

to 

ro 

v v4 

W4 

v\4 

ro 

rsj  c\4 

X <\* 

•r» 

-w  o :< 

ro 

ri 

• 

r u 

>- 

•V4  *^1  *X» 

t>4  ¥ 

4-^ 

v'4 

1 * 

• 

wo 

*t 

IO 

X 

v\l 

X 

- 4 r^ 

~+ 

4-4 

r>- 

• 

J It 

I— 

2r 

c.» 

mikui 

>< 

o ♦♦ 

o »v 

•.  * 

4 

Ui 

a 

K 

1 

UL 

U5  ••• 

» « 

U UI 

i - 

h“ 

1- 

144 

v» 

»»a 

U 

lAJ  VO 

J » 4 Ui 

Vi 

• o 

u> 

• J* 

J 

V.  !• 

u 

a *-d 

m 

4 

% 

• 

ro 

wo 

»15  1 05 

•T. 

4 

04 

|4» 

wo 

r - 

•y,  ‘Xt 

• • 1 

I'l 

Nt 

w • 

*4.' 

1 . 

t.»  I 

n. 

-4  Oj 

1 '‘I 

o 

in 

‘a  r- 

»3J  <X> 

o. 

(V. 

<Xi 

«x.  or.  »x< 

•-  4 

» -4 

» ^ 

•-4 

• •• 

*-4 

— • '•  J 

r i 

• 1 

' J 

• ‘J 

• *J 

v>4  < 1 

• *i 

•i 

1 4 . 

r* . i*“. 

|4"» 

io 

1 4 . 

n . 

IO 

•X.  «• 

'Xi  iV 

vi* 

<X#  05  Vi* 

c.. 

Oi 

•ii  'i* 

Oi 

t5> 

tx* 

til  Ca 

Ca 

•i> 

•X- 

til 

vi* 

tX> 

•Xi 

til 

•ii 

ti» 

•Xt  ‘ii 

•Xi 

a. 

X)  w «v 

E. 


Computational  Analysis 

Using  the  HP-67  standard  curve  fitting  program,  a good 
linear  fit  was  obtained  on  the  four  engine  maximum  range  data. 
There  is  a linear  relationship  between  altitude  and  indicated 
airspeed  for  each  gross  weight  category.  The  coefficient  of 
determination  was  equal  to  1.00  in  all  cases,  indicating  a 
good  fit.  The  following  equations  were  used;  loiter  airspeeds 
are  constant  for  each  category. 

G.W.  (1000  lbs)  Max  range  IAS  Loiter  IAS 


132.5- 

127.5 

y 

= 

270 

— 

x/1000 

220 

127.5- 

122.5 

y 

* 

267 

- 

x/1000 

215 

122.5- 

117.5 

y 

= 

265 

- 

x/1000 

210 

117.5- 

112.5 

y 

= 

262 

- 

x/1000 

205 

112.5- 

107.5 

y 

= 

260 

- 

x/1000 

200 

107.5- 

102.5 

y 

= 

257 

- 

x/1000 

195 

102.5- 

97.5 

y 

= 

255 

- 

x/1000 

190 

97.5- 

92.5 

y 

= 

252 

- 

x/1000 

185 

92.5- 

87.5 

y 

= 

250 

- 

x/1000 

180 

87.5- 

82.5 

y 

= 

247 

- 

x/1000 

175 

1 

in 

• 

CM 

00 

77.5 

y 

= 

245 

- 

x/1000 

170 

77.5- 

72.5 

y 

- 

242 

- 

x/1000 

165 

x = altitude  in  feet  and  y = maximum  range  IAS. 
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USER-CONTROLLED  SIMULATION  OF  APPROACH  AND  LANDING  FOR 


THE  P-3  AIRCRAFT  by  LT  J.  Aiken 

A.  Problem  Statement 

This  program  simulates  an  aircraft  approach  and  landing. 
Specifically,  it  is  a time-step  simulation  of  the  final  five 
miles  of  a precision  approach  for  a Lockheed  P-3  ORION  aircraft. 
The  simulation  is  user-controlled  which  allows  the  user  to  act 
as  pilot  and  make  the  decisions  which  control  the  movement  of 
the  airplane  during  its  final  approach  phase.  The  purpose  of 
the  program  is  to  simulate  accurately  the  flight  of  the  air- 
craft and  display  to  the  operator  his  rate  of  movement  and 
position  resulting  from  his  manipulation  of  the  controls. 


B.  Operation  Analysis 

Relevant  information  on  the  airfield  is  as  follows: 

Runway  8000  ft  (length 

200  ft  (width) 

180  degrees  magnetic  heading 
SEA  LEVEL  elevation 

Approach  TOUCHDOWN  POINT  1000ft  beyond  approach 

threshold 

GLIDE  SLOPE  2.83  degrees 

2 min  18  sec  time  required  at  135  kts 

ground  speed 

FINAL  APPROACH  FIX:  5 miles,  1500ft 

(starting  point) 
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PRICKLING  FjflS  HLAMC 


The  aircraft  weighs  90,000  lbs  with  approach  speeds 
135/121  kts  (approach  flaps/land  flaps) . Note  however  that  no 
provision  is  made  for  changing  the  gear/flap  configuration  so 
it  is  essentially  an  "approach-flap"  landing.  The  simulation 
starts  with  the  aircraft  in  motion:  1500  feet  MSL,  135  kts  IAS, 
650  ft/min  descent  rate,  landing  gear  down  and  approach  flaps. 

The  simulation  allows  the  user  to  select  horsepower  settings, 
nose  attitude,  heading,  wind  direction,  wind  velocity,  and  time 
interval.  At  the  end  of  a time  interval  the  simulation  is 
halted  and  the  critical  flight  information  is  displayed  to  the 
operator,  allowing  him  to  alter  controlling  parameters  and 
continue  the  flight.  The  simulation  continues  in  this  manner 
until  the  aircraft  lands.  Vital  landing  parameters  are  displayed 
and  the  simulation  is  complete.  The  simulation  realistically 
responds  to  control  changes  provided  the  aircraft  is  flown  in 
a somewhat  reasonable  fashion.  Extreme  deviations  and  maneuvers 
other  than  those  required  during  an  approach  are  not  designed 
into  the  program. 
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Computational  Algorithm 

Initialize  the  aircraft  at  the  starting  point. 

Input  time  step,  wind  direction  and  wind  velocity. 

Input  horsepower,  nose  attitude,  heading,  and  number  of 
time  steps  desired. 

Compute  course  deviation. 

Compute  horizontal  acceleration. 

Compute  vertical  acceleration. 

Compute  final  velocity  and  average  vertical  velocity. 
Compute  altitude. 

Compute  final  and  average  horizontal  velocity. 

Compute  distance  remaining  based  on  ground  speed. 

Compute  glide  slope  height  and  deviation  from  glide  slope, 
Check  altitude  less  than  0. 

DSZ  (number  of  time  steps  is  the  counter)  GTO  4 above. 
Display  approach  parameters  after  completing  desired 
time  steps. 

Display  landing  parameters  upon  landing. 

Clear  primary  and  secondary  registers,  GTO  1 for  now 
problem. 


D.  HP-67  Calculator  Program 


1.  User  Instructions 


Instruction 


Enter  program  card 

Clear  primary  and  secondary 
registers 

Initialize 

Enter  time  step 

Enter  wind  direction 

Enter  wind  velocity 

Enter  horsejxjwer 

Enter  nose  attitude 

Enter  heading 

Enter  number  of  time  steps 


Output 


Altitude.  Airspeed  (Packed) 

Descent  rate  (ft/min) 

Above  (+)  or  below  (-) 
glide  slope 

Distance  to  go 
Airborne  (miles) 

Landed  (feet) 


Right  ( + ) or  left  (-)  of 
of  course 


I f altitude  GT  zero  go  to 
Step  4;  make  new  entries 
only  if  change  desired. 

Step  10  must  be  re-entered. 


Input  Key (s) 


se  con  ds 


degrees 


knots 


+ degrees 


degress 


integer 


STO  C 


STO  D 


STO  E 


Output 


30384 


Nose  attitude 


Heading 

flashing 


Alt.  airspeed 


R/S  I Descent  rate 


R/S  Feet  hi/lo 


R/S  + Distance  to 


ft  right/left 


1 


Output 


Input 


ft 

3 

W 


§ 

•H 

4J 

U 


>1 

4J 


2 

•H 

Q 


U 

O 


£ 


T3 

■S 

3 


TJ 

•5 

3 


0) 

I 

a 

<u 

(0 

n 

0 

x 


TJ 

3 

P 

•H 

P 

P 

< 


<u 

to 

o 

z 


cn 

c 

•H 

■o 

(0 

a> 

X 


(0 

ft 

p 

w 


20  210 


20 


10 


& 

o 


w 


a> 

T5 

•H 


2 


o 


« 

p 

4-> 


4J 

3 


T3 

a> 

| 

w 

U 

•H 

< 


£ 


o 


5 

8 


V 

C 

Q) 

O 

W 


0) 

2 


8 

§ 


a 


•P 

W 


45 


<U 

0) 

U 

P 

o 

u 


■P 

«P 


a 


5s 


800 

1.5 

185 

1 

1295 

133 

-570 

-1 

4.4 

63 

183 

1 

1108 

132 

-552 

2.4 

3.  7 

-34 

790 

184 

1 

922 

132 

-560 

6.9 

3.1 

-51 

775 

185 

1 

732 

132 

-580 

6.3 

2.4 

11 

184 

1 

538 

132 

-584 

1.6 

1.8 

- 6 

790 

184 

1 

441 

132 

-575 

0 

1.5 

-14 

184.8 

1 

346 

132 

-574 

-1 

1.2 

9 

795 

2 

184.  3 

1 

252 

131 

-555 

-1 

.85 

13 

700 

3.5 

184 

1 

205 

129 

-559 

-1 

.7 

8 

4 

1 

159 

127 

-552 

-2 

.5 

4 

184.1 

1 

113 

127 

-552 

-2 

.4 

2 

4.3 

184 

1 

67 

127 

-547 

-3 

.24 

- 2 

4.8 

1 

22 

126 

-538 

-3 

.08 

- 7 

6 

1 

4 

126 

-530 

-3 

.02 

- 8 

8 

182.5 

1 

0 

125 

-516 

-2 

40 

-22 

,i 


■ - — 


363$4.0t’  ♦*» 

26.66  i It  w 

216.66  Sill 

. t3t/  v 1 w£ 
800* &C  bSEr. 

1.50  £362 
185.00  6SE2 


1.88 

bOt. 

235. 133 

»♦* 

R.  3 

-570. 66 

**♦ 

R''  s> 

-1.62 

*•* 

R-  S 

4.35 

**♦ 

R/S 

62.92 

**» 

183.08  £622 

i . 88  OSSv 
1106.132  *<M 


■552.06 

R/S 

2.42 

R/S 

*♦* 

3.71 

R/S 

$** 

-34.34 

R/S 

**» 

790.60 

eSfcn 

164.88  £362 

1 . 88  £S6u 

» 22.132 

*** 

560.46 

R/S 

**♦ 

-51.46' 

"5.66 

135.88 

1.00 

732.132 


R-S 

1.00  kSEZ 

-555.34  **? 

67.127  **♦ 

R/S 

R/S 

-1.25  **« 

-548.74  *44 

R-'S 

R-  S 

0.85  *»*  ■ 

-2. 72  44* 

R-S 

R-S 

12.69  *** 

0.24  *+4 

5.k>0  STCC 

R-  2 

700.00  i-iin 

-2.43  *** 

5.50  bSe-S 

4mSd  655* 

164.66  bS£J 

1.00  bSEZ 

1.00  fssr 

22.128  44* 

205.129  *** 

R/S 

R'  S 

-5 36.29  44* 

-559.51  *** 

R- 3 

R/S 

-2.76  *4* 

-1.27  44* 

R-  £ 

R-S 

0.08  444 

0.69  *** 

R'  c 

R-S 

-6.71  #*< 

6.41  **4 

2.00  ST2Z 

4.00  £S£E 

6.06  CSEt 

1 . 09  SSf'Z 

1.00  £SEZ 

159.127  «** 

4.126  4*4 

R/£ 

R-S 

-552.84  *** 

-530.41  44* 

R-S 

R-S 

-1.65  *** 

-2.61  44* 

R-  i 

R-S 

0.54  *** 

0.62  444 

R-S 

R.  S 

4.13  44* 

-8.42  4*4 

104.10  bS£2 

8. 62  65EE 

1.00  bSEZ 

182.50  8S£C  ' 

113.127  *+•»■ 

1.00  6SEZ 

R-S 

0.125  44* 

-551.51  #44 

R.  S 

R-S 

-516.64  #*•» 

^ ^ 

R/S 

R/S 

-2.24  #** 

0. 59  *** 

R-S 

R-S 

46.50  4*4 

1.85  *4* 

R-S 

4.36  CS6t 

-22.16  4*4 

184.00  £582 

3.  Program  Storage  Allocations  and  Program  Listing 
Registers : 


RO  : 

Right/left  of  course 

81: 

Horti zonal  force 

Rl: 

Horsepower 

S2  : 

Vertical  force 

R2 : 

Nose  attitude 

R3: 

Horizontal  velocity 

R4  : 

Vertical  velocity 

R5 : 

Altitude 

R6 : 

Horizontal  acceleration 

R7: 

Vertical  acceleration 

R8 : 

Distance  remaining 

R9 : 

Glide  Slope  Altitude 

RA: 

DELH 

RB: 

Heading 

RC: 

Time  step 

RD: 

Wind  direction 

RE: 

Wind  velocity 

— 

RI: 

Number  of  time  steps 

Initial  Flag  Status; 


0: 

OFF, 

Unused 

2: 

OFF, 

Set  ON  upon  landing 

1: 

OFF, 

Unused 

3: 

OFF, 

Unused 

User  Control  Keys: 

A:  Horsepower  a: 
B:  Nose  attitude  b: 
C:  Heading  c: 
D:  Time  steps,  start  computation  d: 


E: 


e:  Initialize 


001 

*L?Lh 

21  11 

002 

$731 

35  01 

eoj 

R $ 

51 

004 

*LEL6 

21  12 

005 

$ to: 

35  v2 

m 

R'S 

51 

997 

tLBLC 

21  13 

005 

sros 

35  12 

00? 

R-S 

51 

010 

»LELe 

21  1?  15 

011 

3 

4b 

02 

012 

<J 

02 

01  J 

? 

02 

014 

STP3 

35  03 

015 

1 

01 

010 

1 

01 

01? 

CHS 

^ "l 

010 

ST04 

35  04 

01? 

1 

01 

020 

5 

C5 

021 

0 

00 

022 

0 

00 

023 

S7P5 

35  05 

024 

3 

03 

025 

0 

DC 

020 

3 

03 

927 

0 

0? 

92S 

4 

04 

929 

$70? 

35  02 

030 

R $ 

51 

031 

*LFLP 

21  14 

032 

STOJ 

35  4? 

033 

*Lf  L? 

21  0? 

034 

RCL6 

36  12 

03  5 

1 

01 

030 

0 

OS 

07? 

0 

9C 

07? 

- 

-45 

03? 

SIN 

41 

040 

J 

4> 

02 

041 

3 

03 

042 

0 

00 

043 

X 

-35 

044 

I 

01 

045 

0 

0? 

04t> 

0 

00 

04? 

RClD 

3?  14 

04? 

- 

-45 

04? 

s;n 

41 

050 

RCLE 

Jo  15 

051 

1 

01 

• 052 

• 

-02 

053 

6 

00 

054 

9 

0? 

055 

x 

-35 

05? 

V 

-35 

05? 

t 

-55 

05? 

WLC 

30  17 

05? 

A 

-35 

960 

ST*fl 

35-55  CO 

INITIALIZE 


SPEEDS 


ALTITUDE 


DISTANCE 


COMPUTE 


DR I IT  RATES 


COURSE 


DEVIATION 


001 

CSFi 

23  01 

002 

ST  00 

35  00 

007 

CSF2 

23  02 

004 

st  or 

35  07 

005 

r:l4 

30  04 

000 

RCL  7 

30  0? 

007 

RCLC 

3o  13 

06? 

X 

- -35 

00? 

S7*4 

35-55  04 

070 

\;v 

-41 

071 

K CL  4 

30  04 

072 

♦ 

-55 

673 

2 

02 

074 

T 

s « 

“4l** 

075 

RCLC 

5t>  i 3 

070 

X 

-35 

077 

ST+5 

35-55  05 

07? 

RC-3 

30  02 

07? 

RCLb 

30  06 

0SO 

RCLC 

30  13 

0?J 

X 

-35 

982 

$7+3 

35-55  03 

983 

xrv 

-4.' 

9S4 

RCL3 

30  07 

0?5 

♦ 

-55 

0?0 

2 

02 

98? 

* 

-24 

989 

RCLO 

36  14 

9S9 

1 

01 

990 

? 

02 

991 

0 

00 

992 

- 

-45 

997 

HPS 

10  31 

994 

COS 

42 

995 

RCLE 

30  15 

996 

1 

01 

997 

s 

-02 

998 

6 

00 

999 

0 

0? 

100 

A 

-35 

101 

\ 

-35 

102 

- 

-45 

103 

RCLC 

36  17 

104 

\ 

-35 

105 

CHS 

-22 

106 

RCL? 

36  02 

107 

♦ 

-55 

10? 

$70? 

35  0? 

10? 

a 

-02 

110 

0 

00 

111 

4 

04 

112 

.4 

0? 

1 13 

A 

-35 

114 

?T.T? 

35  03 

115 

CHS 

110 

RCL  5 

36  05 

11? 

♦ 

-55 

11? 

ST  Oh 

35  11 

11? 

RCL  5 

30  05 

120 

.V  0" 

lb-41 

121 

CTO? 

22  0? 

122 

EC2I 

Jo  25  4 L 

COMPUTE 

HORIZONTAL  ACCEL 


COMPUTE 

VERTICAL  ACCEL 


COMPUTE 

VERTICAL 

VELOCITY 


COMl’UTE 


ALTITUDE 


COMPUTE 

HORIZONTAL 


VELOCITY 


COMPUTE 
C.LI  DE-SLOPE 


HE  1 C.HT 


DEVIATION 


OLI DE-SLOPE 


CHECK  ALT 


UXIP  CONTROL 


124 

fid? 

21  87 

169 

iLBcl 

cl  8 . 

125 

fC  15 

38  85 

178 

f:l2 

38  82 

128 

m 

18  .'4 

171 

-82 

127 

12$ 

FCL3 

1 

38  83 

6. 

ALTITUDE 

172 

17? 

4 

; NS 

84 

"CC 

COMPUTE 

129 

iso 

6' 

0 

w 

88 

8c 

DECIMAL 

174 

175 

> 

EE.V 

_ ■ c 

"C  c 

HORIZONTAL 

131 

8 

8t 

AIRSPEED 

1 78 

3 

83 

132 

13? 

4 

f 

-c4 

-55 

(PACKED) 

177 

17? 

CHS 
f CL  I 

-22 
38  81 

FORCE 

134 

Pi'c3 

-83  83 

17? 

K 

-35 

135 

R S 

51 

188 

♦ 

_ r r 

w «• 

136 

PSP2 

-83  82 

1?1 

v 

-62 

13." 

P2L4 

38  04 

RATE 

182 

4 

Of 

138 

6 

88 

OF 

183 

- 

-45 

13? 

e 

88 

184 

FiS 

18-5. 

148 

A 

-35 

DESCENT 

185 

Fdl 

38  81 

141 

R S 

51 

1S6 

v * V 

-41 

142 

143 

RCcR 

R $ 

38  11 

51 

HIGH/LOW 

IS 7 
183 

ST01 

,\:v 

35  81 

-41 

COMPUTE 

144 

145 

RCL8 

F2A 

38  8c 

18  23  82 

18? 

1?8 

f;s 

-45 

18-51 

HORIZONTAL 

146 

147 

148 

6TQ5 

6 

8 

22  85 

88 

88 

DISTANCE  TO  GO 

(FEET  IF 

1?1 

192 

197 

RCL8 

5 

38  88 
85 
-24 

ACCELERATION 

149 

158 

7 

I 

87 

LANDED) 

194 

195 

FTH 

-55 

24 

151 

152 

*L8:f 

-24 

21  85 

196 
19  7 

tLBL  2 

21  82 
-82 

153 

R'S 

51 

198 

0 

00 

154 

155 

PC  id 
R-S 

38  88 

51 

LEFT/ RIGHT 

199 

208 

0 

1 

80 

81 

COMPUTE 

156 

tLBl? 

21  8? 

201 

RCL1 

38  81 

157 

158 
15? 
1(6 

RCL? 

RCL5 

CHS 

* 

38  8t 

38  85 

-££ 

82 

LANDED.  . . 

ADJUST  DISTANCE 

282 

283 

284 
205 

\ 

/ ■ 

8 

-35 

-82 

0? 

-45 

VERTICAL 

FORCE 

161 

162 

8 

\ 

88 

-35 

TO  ZERO  ALT. 

288 

287 

R2w2 

38  82 
* 82 

167 

1(4 

4 

sro8 

-55 

35  88 

20S 

209 

8 

5 

88 

05 

165 

8 

88 

218 

,\ 

-35 

166 

ST  05 

35  85 

SET  ALT  - 0 

211 

* 

-55 

16? 

ST2 
iT  07 

18  21  82 

SET  FLAG  2 

212 

p:s 

18-51 

1(8 

22  87 

213 

r:l: 

36  82 

214 

215 
218 

V-V 

-41 

STC2 

x:v 

35  02 
-41 

COMPUTE 

217 

218 

F7S 

-45 

J8-51 

VERTICAL 

213 

228 

RCL  7 

5 

38  07 
85 

ACCELERATION 

221 

4> 

-2  4 

222 

♦ 

r r 

223 

PTH 

24 

224 

R.  $ 

5i 

50 


Mathematical  Analysis 


The  variable  names  used  are: 

VX  Horizontal  Velocity 

VY  Vertical  Velocity 

VW  Wind  Velocity 

DRA  Drift  Rate  Due  to  A/C  HDG 

DRW  Drift  Rate  Due  to  Wind 

TDR  Total  Drift  Rate 

GS  Ground  Speed 

AX  Horizontal  Velocity 

AY  Vertical  Velocity 

D Distance  to  go 

H Glide  Slope  Altitude 

DELH  Deviation  from  Glide  Slope 

ALT  Altitude 

t Time  Step 


51 


LANDING  DISTANCE  ADJUSTMENT 


UPON  LANDING 


J 

■ 

I 

■ 

■ 


Two 

situations  are 

illustrated 

above . 

1. 

Aircraft 

above 

glide 

slope 

(lands 

long) . 

2. 

Aircraft 

below 

glide 

slope 

( lands 

short)  . 

An  adjustment  is  required  because  the  aircraft  is  allowed  to 
descend  until  the  end  of  the  time  step.  On  the  landing  time  j" 

step  the  aircraft  will  have  descended  below  zero  altitude.  | 

The  adjustment  in  either  case  is:  adjustment  = H/(tan  0). 

II 


FORCES,  ACCELERATIONS 


FX 

= - e + 

. 001HP  - 

.4 

FY 

= . 001HP  - .8  + 

.050 

AX' 

= (FX1  - 

■ FX° ) + 

. 2AX° 

AY' 

= (FY1  - 

+ 

o 

X 

Cm 

. 2AY° 

VX1 

II 

< 

X 

o 

+ 

AX  ( t) 

VYf 

+ 

o 

£ 

n 

AX  ( t) 

AND  VELOCITIES 

FX  = horizontal  force 

FY  = vertical  force 

0 = nose  attitude 

Superscripts  denote  different 
time  steps 

Subscripts  denote  start  and  end 
of  a time  step. 


Although  the  author  is  not  an  aeronautical  engineer,  it 
was  felt  that  his  understanding  of  the  basic  laws  of  physics  com- 
plemented by  considerable  pilot  experience  would  serve  sufficiently 
to  accomplish  the  goals  of  this  project.  The  formulae  for  force 
and  acceleration  were  arrived  at  after  testing  several  trial 
formulae  on  experienced  P-3  pilots.  The  unanimous  opinion  was 
that  the  current  program  enables  the  simulation  to  closely  model 
the  actual  flight  characteristics  of  the  aircraft. 
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VI.  FLIGHT  CREW  MANAGEMENT  USING  THE  HP-97  by  LT  Kenneth  W.  Peters 


A.  Problem  Statement 

A flight  crew's  most  recent  landing  day  and  time  is  known. 
Using  requirements  for  crew  rest  and  postflight  and  preflight 
duration,  compute  when  the  flight  crew  will  be  available  for  take- 
off again.  For  planning  and  scheduling  purposes,  list  crews  in 
order  of  availability.  For  required  onstation  times  compute 
takeoff,  onstation,  offstation  and  landing  times  for  a given 
number  of  flights.  Determine  if  flight  crews  will  be  available 
to  meet  this  schedule. 


B.  Operational  Analysis 

When  planning  an  operation  requiring  scheduling  of 
several  flight  crews,  crew  availability  must  be  considered. 
Accurate  and  easily  understood  crew  records  are  necessary  to 
meet  both  operational  and  safety  requirements. 


C.  Computational  Algorithm 
1.  Flight  crew  availability 

a.  Enter  required  postflight  to  preflight  crew  rest  time. 

b.  Enter  crew  number  and  their  most  recent  landing  day  and  time. 

c.  Compute  the  crew's  earliest  possible  takeoff  day  and  time 
using  one  hour  for  postflight  and  three  hours  for  preflight. 
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PRECEDING  PAGE  bLAlflC 


List  crews  in  order  of  availability 

a.  Enter  number  of  flight  crews  = N. 

b.  Compare  crew  N availability  with  crew  (N-l)  availability. 
Store  crew  number  and  availabilicy  of  crew  which  can  take- 
off soonest.  Compare  this  crew  with  crew  (N-x) , 

(x  = 1 , 2 , 3, 4 , . . . , (N-l) ) , and  store  the  crew  which  is 
available  the  soonest. 

c.  When  most  available  crew  has  been  determined,  add  10,000 
days  to  its  takeoff  availability  date  and  increment 
counter. 

d.  Repeat  Steps  b and  c until  all  crews  have  been  listed, 
i.e.  counter  = N. 

e.  Restore  crew  availability  data  by  subtracting  10,000  days 
from  each  crew's  availability  date. 

Operational  schedule. 

a.  Enter  required  number  of  flights  from  a particular  base 
and  the  gap  (+)  desired  for  onstation  coverage. 

b.  Compute  takeoff  day  and  time,  onstation  day  and  time, 
off station  day  and  time,  and  landing  day  and  time. 

c.  Check  number  of  flights  versus  counter.  Repeat  b as 
required. 

d.  Compare  print-out  of  required  takeoff  times  with  crew 
availability  listing  (see  Part  2 above) . 
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D.  HP-67/97  Calculator  Program 
1.  User  Instructions 


Step 

Instruction 

Input 

Key (s) 

Output 

1. 

Enter  program  card 

2. 

Compute  flight  crew  availability 

a. 

Enter  crew  rest  time 

HH.MM 

fA 

b. 

Enter  crew  rest  data: 

Crew  # 

1 < # £ 19 

+ 

Postflight  time 

HH.MM 

+ 

Preflight  time 

HH.MM 

t 

Landing  Julian  day  and  time 

DAY.HHMM 

A 

Crew  # 

Availability 

DAY.HHM 

NOTE:  If  change  in  year  will 

occur,  Julian  date  can 
be  entered  as: 

YRDAY.HHMM 

3. 

List  crews  in  order  of  avail- 

Crew # 

ability. 

Availability 

Enter  # of  crews  = N 

1 < N < I9 

B 

DAY.HHMM  (re- 
peated N times) 

4. 

Compute  operational  flight 

schedule 

a. 

Enter  # of  flights  and  gap 
in  on station  coverage: 

# of  flights  = N 

1 < N < 25 

+ 

gap  in  coverage 

HH.MM 

fC 

b. 

Compute  schedule 

onstation  day  and  time 

DAY.HHMM 

t 

one-way  transit  time 

HH.MM 

+ 

mission  time 

HH.MM 

c 

Flight  # 

Takeoff  DAY.HHMM 
Onsta  DAY.HHMM 
Offsta  DAY.HHMM 
Land  DAY.HHMM 

5. 

(Optional)  Store  crew 
availability  data  on 
magnetic  card: 

Set  W/P RGM- RUN  switch 
to  RUN 

f W/DATA 

To  reload  data: 

Set  W/P RGM- RUN  switch 
to  RUN  and  load  data  card 
obtained  above 

1 

2 . Sample  Problem 


a.  Input  the  following  data: 
Crew  rest  = 15  hours 


Crew  # 
1 
2 

3 

4 


Postflight 

1 

1 

1 

1 


Preflight 

3 

3 

3 

3 


b. 


Compute  crew  availability. 
Answer : 


Crew  # 
1 
2 

3 

4 


Availability 

3.1830 

2.0700 

3.1300 

2.1032 


c.  Input  the  following  crew  status: 


Landing 

2.2330 

1.1200 

2.1800 

1.1532 


Crew  # (in)  Register  # Availability 


1 

1 

3.1830 

2 

2 

2.0700 

3 

3 

3.1300 

4 

4 

2.10  32 

5 

5 

2.0958 

6 

6 

2.2020 

7 

7 

3.1930 

8 

8 

5.0730 

9 

9 

4.0930 

10 

10 

6.0430 

11 

11 

4.0345 

12 

12 

2.0730 

13 

13 

8.0130 

14 

14 

5.0315 

15 

15 

3.0545 

16 

16 

3.1300 

17 

17 

9.0955 

18 

18 

2.1030 

60 


d.  Output  listing  of  crews  in  order  of  availability. 


Answer:  2,  12,  5,  18,  4,  6,  15,  3,  16,  1,  7,  11,  9, 

14,  8,  10,  13,  17. 


e.  Develop  an  operational  flight  schedule  consisting  of 
six  (6)  flights  with  a zero  (0)  gap  in  onstation 
coverage.  The  first  onstation  time  is  2.1830,  it 
takes  one  (1)  hour  and  fifty  (50)  minutes  for  a 
one-way  transit,  and  a flight  from  this  base  has 
nine  (9)  hours  of  total  mission  time. 

Answer: 


Flight  # 

Takeoff 

Onsta 

Of fsta 

Land 

1 

2.1640 

2.1830 

2.2350 

3.0140 

2 

2.2200 

2.2350 

3.0510 

3.0700 

3 

3.0320 

3.0510 

3.1030 

3.1220 

4 

3.0840 

3.1030 

3.1550 

3.1740 

5 

3.1400 

3.1550 

3.2110 

3.2300 

6 

3.1920 

3.2110 

4.0230 

4.0420 

f.  Develop  an  operational  flight  schedule  consisting 

of  two  (2)  flights  with  a two  (2)  hour  gap  in  onstation 
coverage.  The  first  onstation  time  is  2.1830,  with 
the  same  transit  and  mission  time  as  in  (e)  . 
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Answer : 


Flight  # Takeoff  Onsta  Offsta  Land 


2.1640  2.1830  2.2350  3.0140 


3.0000 


3.0150  3.0710  3.0900 


g.  Develop  an  operational  flight  schedule  with  the  same 


parameters  as  in  (f) , except  with  a negative  two  (-2) 


hour  gap  in  onstation  coverage. 


Answer : 


Flight  # 


Takeoff 


Onsta 


Offsta 


Land 


2.1640 


2.1830  2.2350 


3.0140 


2.2100 


2.2250  3.0410 


3.0600 


h.  With  one  flight  from  Base  A with  the  first  onstation 


time  as  above  and  one  flight  from  Base  B to  the  same 


operational  area,  compute  the  flight  schedule.  One- 


way transit  time  from  Base  B is  two  (2)  hours,  but 


mission  time  is  now  ten  (10)  hours.  Zero  (0)  gap 


in  coverage  is  desired, 


Answer : 


Flight  # 


Takeoff  Onsta 


Offsta 


Land 


2.1640  2.1830  2.2350 


3.0140 


2.2120  2.2350  3.0450  3.0720 
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With  four  (4)  flight  crews  with  crew  availability 
as  follows: 


Crew  # 


Availability 


1 

2 

3 

4 


3.1830 

2.0700 

3.1300 

3.1700 


determine  if  they  are  available  to  meet  an  oper- 
ational flight  schedule  with  the  following  con- 
ditions : 

4 flights 

2 hour  gap  in  onstation  coverage 
3.0100  is  the  first  onstation  time 
1 hour  and  30  minute  one-way  transit 
10  hour  mission  time 

Answer:  No,  a crew  is  not  available  to  meet  the 

second  scheduled  takeoff  time. 
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Examples  a,  b 


4 5.  Ol'CO  I'.'tv, 

i # oofi'  £h  r .* 

1.0000  EN71 

3.0000  EM  f 
2.2330  Oocn 

1.  »•« 
3.1330  «M 

2.0000  ENT* 

1.0000  EM? 
J.OPO0  EM? 

1.1200  i$Sn 

2.  m» 
2.0700  M» 


J • 

0000 

EM; 

1. 

0000 

EM? 

3. 

0000 

ENT  r 

2. 

1300 

03E* 

3. 

hi 

3. 

J 300 

Ml 

4. 

0000 

ENT  * 

1. 

0000 

ENT 

3. 

0000 

EM* 

1. 

1532 

feibr 

4. 

<M 

2. 

1032 

IM 

Examples  c,  d 


IS. 00 00 

03SE 

15. 

M* 

•» 

Ml 

3.0545 

Ml 

2.0700 

IM 

J • 

Ml 

12. 

Ml 

3.1300 

Ml 

2.0730 

Ml 

lt>. 

Ml 

5. 

Ml 

3. 1 300 

Mi 

2.095S 

Ml 

1. 

IM 

13. 

Ml 

3. 1 330 

1 1 1 

2.1030 

Ml 

7. 

Ml 

4. 

Ml 

3.1930 

Mi 

2. 1032 

Ml 

11. 

• M 

0. 

Ml 

4.0345 

• M 

2.2020 

• M 

S.  «M 

4.0930  *»» 


14.  mi 
'.0315  mi 


b.  mi 
5.0730  mi 


10.  mi 
3. 0430  mi 


13.  m» 
S.0130 


1 1 . *»  i 

.>5.0*55  in 


Example  e 


Example  f 


6.  0000 

El > : * 

O.OtYt 

3. I030 

e,n*: 

1 . 5006 

EMf 

?. 0600 

oseo 

1 . oooooooo 

1*4(5 

»»* 

bi  • 1 650 

3.3350 

■**» 

3.0140 

»** 

2.  oooooooo 

*»l 

2.2200 

**t 

2.2350 

*»* 

3.0510 

*** 

3.0700 

*»* 

3. 00000000 

3.6320 

1*1 

3. 0510 

*»». 

3. 1 030 

t*» 

3. 1220 

»»* 

4.00000000 

3. 0640 

*♦» 

3. 1030 

m 

3.1550 

«#  1 

3. 1 740 

»»* 

5. 60000000 

**» 

3. 1 400 

tit 

3.1550 

tit 

3.3110 

ttt 

3. 2300 

tit 

t. 00 060000 

lit 

3. 1 926 

tit 

3.2110 

in 

4.0230 

tti 

4.0420 

Ml 

2.vv06  Eh i ' 
0 . 0000  Oof. 

2. 1 f 30  th~  ‘ 
1.5000  ENT: 
0.0000  0SEJ 
1 . 00660606  » * * 
£ - 1 1>  4 O »*» 
3 . 1 636  * * * 
3.3350  »»i 


3. 00000000 

Ml 

3. 0000 

*»* 

3.0150 

Ml 

3.6710 

lit 

3. 0906 

li » 

Example  c 


3.0000 

EHT! 

-1.0000 

Oof; 

3.1330 

Eh:: 

1 . 5000 

est; 

?. 0000 

Off  J 

1 . 00000060 

m 

2. 1646 

»**' 

2. 1S30 

9M 

2. 3350 

3.0140 

3. 00000600 

Id 

2.2106 

»♦» 

2.2256 

it. 

3. 0410 

»n 

3.  Ot'OO 

im 

Example  h 


4.6C00 

6666 

£ . 

444 

2.0700 

*♦  ? 

3. 

(U 

' 

3.  1366 

44* 

4. 

4*4 

3.17(90 

)M 

1. 

IM 

3.1638 

444 

4.0000  En7; 

2.0000  666: 

3.0100  EHT! 

1.3000  Ef 1 7 1 

10.0000 

bit  J 

1 . 00000000 

*1  f 

2.2336 

M4 

3.0100 

• 44 

3.0606 

44* 

3.0630 

44* 

2.  00006600 

**  « 

3. 0630 

*4* 

3.1O00 

444 

3. 1760 

44* 

3.1830 

444 

3. 00000068 

3. 1 730 

44* 

3.1f>00 

44* 

4.0266 

4.6330 

44* 

4.60006000 

44* 

4.0230 

444 

4.0406 

44* 

4.1100 

44* 

4.1230 

441 
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3.  Program  Storage  Allocation  and  Listing 


Registers : 


RO 

Crew 

rest 

SO 

Crew 

10 

availability 

Rl 

Crew 

1 

availability 

SI 

Crew 

11 

availability 

R2 

Crew 

2 

availability 

S2 

Crew 

12 

availability 

R3 

Crew 

3 

availability 

S3 

Crew 

13 

availability 

R4 

Crew 

4 

availability 

S4 

Crew 

14 

availability 

R5 

Crew 

5 

availability 

S5 

Crew 

15 

availability 

R6 

Crew 

6 

availability 

S6 

Crew 

16 

availability 

R7 

Crew 

7 

availability 

S7 

Crew 

17 

availability 

R8 

Crew 

8 

availability 

S8 

Crew 

18 

availability 

R9 

Crew 

9 

availability 

S9 

Crew 

19 

availability 

RA 

Landing  data;  # of  crews 

RB 

First 

onsta  time;  crew  #;  takeoff 

time 

RC 

Mission  time;  flight 

counter 

RD 

Gap  between  onstation 

periods;  ith 

availability 

RE 

Flight 

counter;  crew 

counter;  one-way 

transit 

0; 

Unused 

2: 

OFF,  day  correction 

1: 

Unused 

d : 

OFF,  hour  correction 

User 

Controlled 

Keys 

A; 

Crew  # t , 

postflight  t,  preflight  t,  land  =»  compute  availabi 

B: 

crew  # f ; 

compute  listing 

C: 

onsta  t,  one-way  transit  t. 

mission  time  =»  flight  schedule 

D: 

unused 

E: 

unused 

a: 

crew  rest 

d:  unused 

b: 

unused 

e : unused 

c: 

flight  t , 

gap 
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" ■ 11  ' 1 ■ 

eei 

tLFuc 

21  16  12 

Input  and  Store  # of  flights  and  desired  gap 

e*2 

st  or 

35  14 

083 

R4 

-31 

ST  08 

35  11 

005 

RTH 

24 

m 

*LFLC 

21  13 

ee? 

CF2 

16  22  82 

Compute  flight  schedule 

80S 

CF3 

16  22  03 

Input  and  store  onstation  day  and  time 

809 

STCC 

35  13 

Input  and  store  one-way  transit  time 

818 

R4 

-31 

Input  and  store  total  mission  time 

81 1 

STOE 

35  15 

812 

Ri 

-31 

813 

ST  OB 

35  12 

814 

GSB8 

23  60 

81 5 

RCLE 

36  15 

Unpack  ddd.hhmm  to  hh.mm 

816 

CHS 

“22 

81 7 

HHS* 

16-55 

Compute  takeoff  time 

010 

X .0'’ 

16-45 

Check  for  time  greater  than  24  hours 

019 

GSEe 

23  16  If 

828 

EEs 

-23 

021 

2 

82 

822 

~ 

“24 

823 

RCLE 

36  12 

824 

IN7 

16  34 

025 

F3r 

16  23  83 

826 

GSBel 

23  16  14 

827 

•f 

-55 

Takeoff  ddd.hhmm 

828 

STOE 

35  12 

829 

1 

01 

838 

*LBL9 

21  89 

831 

832 

STCJ 

DBFS 

35  46 
-63  8S 

Start  loop 

833 

FFTX 

-14 

Print  flight  # 

834 

RCLB 

36  12 

835 

f 

16  23  62 

836 

CSE2 

23  62 

837 

0SF4 

-63  84 

838 

FRTX 

-14 

Print  takeoff  ddd.hhmm 

839 

STOE 

35  12 

848 

ST08 

35  66 

841 

GSP0 

23  66 

842 

RCLE 

36  15 

843 

HHS* 

16-55 

Compute  and  print  onstation  ddd.hhmm 

844 

GSEO 

23  14 

845 

CSBE 

23  15 

846 

GSB0 

23  66 

847 

RCLE 

36  If 

846 

RClE 

36  If 

849 

hHS* 

16-55 

050 

CHS 

“22 

051 

RCLC 

36  13 

05  2 

HHS* 

16-55 

853 

HHS* 

16-55 

Compute  and  print  offstation  ddd.hhmm 

854 

csbo 

23  14 

055 

CSBE 

23  15 

056 

6SB8 

" 23  68 

05? 

RCLE 

36  If 

856 

HI IS* 

16-55 

Compute  and  print  landing  ddd.hhmm 

05? 

QSBD 

23  14 

060 

CSBE 

23  15 

061 

SFC 

16-11 

68 


Compute  next  takeoff  ddd.hhmm 


062 

PC  LB 

36  12 

863 

6388 

23  6C 

86* 

f:l: 

36  i 3 

865 

hns* 

16-55 

866 

r :le 

36  15 

0t>' 

RCLE 

3b  i5 

868 

HNS* 

16-55 

869 

CHS 

8~e 

HNS ♦ 

16-55 

871 

PC  LD 

36  14 

872 

HI IS* 

16-55 

873 

GSBD 

23  14 

874 

RCLB 

36  12 

875 

IHT 

16  34 

876 

* 

-55 

877 

STOB 

35  12 

876 

RCLH 

36  11 

879 

1 

81 

8S8 

* 

-55 

88 1 

J 

81 

882 

RCL1 

36  4 S 

863 

* 

-55 

884 

v=v? 

16-33 

885 

R s 

51 

886 

GTC9 

22  89 

887 

*LEL& 

21  14 

888 

62 

689 

4 

84 

898 

XtT' 

16-35 

891 

GSBS 

23  86 

092 

x:y 

-41 

893 

EEX 

-27 

894 

2 

82 

09  5 

T 

-24 

896 

RTH 

24 

897 

*LBLS 

21  06 

898 

CHS 

-22 

899 

HNS + 

16-55 

188 

SF2 

16  21  82 

181 

xrv 

-41 

1 82 

RTH 

24 

183 

*LBLE 

21  15 

184 

RCL8 

36  66 

105 

IHT 

16  34 

186 

F2 ? 

16  23  82 

187 

CSB2 

23  62 

186 

♦ 

-55 

109 

PR  TS 

-14 

118 

STOB 

35  66 

lit 

RTH 

24 

112 

»L  BL2 

21  82 

113 

1 

61 

114 

♦ 

-55 

115 

RTH 

24 

11 6 

*LBL0 

21  08 

117 

FPC 

16  44 

118 

EES 

119 

a 

c: 

1 20 

X 

-35 

Increase  counter 


Check  for  exit  from  loop 


Subroutine  to  correct  for  time  greater 
than  24  hours 


Subroutine  to  correct  time 


Output  subroutine 

Print  onstation,  offstation,  and  land 


Subroutine  to  correct  date 


Subroutine  to  change  .hhmm  to  hh.mm 


121 

Pwh 

it 

122 

$LPLe 

21  16  15 

12? 

*'  ? 

02 

124 

4 

64 

125 

Hlr* 

16-55 

126 

rfo 

16  21  6? 

12? 

PTH 

24 

123 

*LE^<J 

21  16  14 

129 

1 

a: 

136 

- 

-45 

131 

RTN 

24 

1 32 

*LFU 

21  16  11 

1 33 

sroo 

35  60 

134 

CF2 

16  22  02 

13  5 

RTN 

24 

136 

$LBLh 

21  11 

13? 

stop 

35  11 

138 

Ft 

16-31 

139 

STOI 

35  46 

148 

F4 

-31 

141 

FRC 

16  44 

142 

£EX 

-23 

143 

4. 
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144 

X 

-35 

145 

HNS* 

16-55 

146 

HNS* 

16-55 

14 ? 

FCLB 

36  02 

148 

HNS * 

16-55 

149 

3SF0 

23  14 

158 

ST01 

35  45 

151 

RCLH 

36  11 

152 

INT 

16  34 

153 

F2? 

16  23  02 

154 

6SB2 

23  02 

155 

ST*i 

35-55  45 

156 

RCL1 

36  46 

157 

DSF8 

-63  00 

158 

FPTX 

-14 

15? 

FCLi 

36  45 

160 

CSP4 

-63  04 

161 

PFry 

-14 

162 

SPC 

16-11 

163 

RTN 

24 

Subroutine  to  correct  for  negative  time 

Subroutine  to  correct  date 

Store  crew  post  flight  to  pre flight  rest  time 

Compute  crew's  earliest  possible  takeoff 
Store  latest  landings  ddd.hhiran 
Store  crew  # 

Compute  crew  I's  takeoff  availability 

Print  crew  # 

Print  crew's  ddd.hhmm  takeoff  availability 


List  flight  crews  in  order  of  availability 
Store  # of  crews  = N 
Zero  counter 


Begin  loop 


Recall  crew  # N availability 


Recall  crew  # N-l  availability 


Compare  availability 


Store  earliest  ddd.hhmm 


Swap  registers 


Output  routine 
Print  crow  tt 


Print  crew  availability  ddd.hhmm 


Add  1000  to  day--i.e.  create  an  artificially 
large  date 


Check  for  end  of  loop 


Subtract  10,000  from  date--i.e.  correct 
artificial ly 


VII. 


TARGET  MOTION  ANALYSIS  (TMA)  OF  A BEARINGS-ONLY  TARGET 


FROM  A MOVING  PLATFORM  by  LT  P.  W.  Marzluff  and 
LT  R.  C.  Pilcher 
A.  Problem  Statement 

Bearings  to  a target  either  stationary  or  moving  with 
constant  course  and  speed  are  available  from  a non-stationary 
tracking  platform.  Determine  the  target's  range,  course  and 
speed. 


B.  Operational  Analysis 

The  four  bearing  TMA  technique  used  in  this  program  requires 
a minimum  of  four  target  bearing  observations  taken  during  a 
minimum  of  two  tracking  legs.  A target  bearing  observation 
must  be  made  and  entered  for  the  time  corresponding  to  the 
initiation  of  own  ship  course  or  speed  change.  Exact  target 
bearing  observations  of  090°  and  270°  require  the  addition  of 
0.1°  to  the  observed  value  to  avoid  infinite  computational 
values.  When  the  tracking  problem  carries  into  a new  day,  the 
previous  day's  time  scale  must  be  continued  (i.e.  a time  of  0010 
on  the  second  day  must  be  enetred  as  2410)  . 

Own  ship  and  target  information  is  entered  on  card  1. 
Estimation  of  target  parameters  begins  on  card  1 and  is  completed 
on  card  2.  Entering  supplemental  target  or  own  ship  informa- 
tion and  generating  a new  estimate  again  requires  the  use  of 
card  1 and  then  card  2 . 

The  accuracy  of  the  estimates  are  dependent  on  the  accuracy 
of  the  inputs,  principally  the  target  bearing,  the  magnitude 
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of  course  or  speed  change,  and  the  number  of  observations  made. 
For  a more  complete  examination  of  the  character  of  the  esti- 
mates, see  Reference  1. 


C.  Computational  Algorithm 

1.  Input  own  ship's  course  and  speed.  Calculate  and  store 
velocity  components. 

2.  Input  time,  t^,  and  observed  target  bearing,  . 

Calculate  elasped  time  since  the  first  observation, 

At^,  and  tan  . Calculate  and  store  the  matrix  values. 

3.  When  own  ship  changes  course  or  speed  enter  t^  and  B^ 
observed  at  the  time  the  course  or  speed  change  was  made. 

Enter  new  own  ship  course  and  speed  prior  to  entering 
the  next  bearing  observation. 

4.  When  at  least  four  target  bearing  observations  have  been 
entered  (bearings  taken  on  a minimum  of  two  tracking  legs) , 
estimation  of  target  range,  course,  and  speed  can  be  made. 

5.  Additional  target  and  own  ship  information  can  be 
entered  and  new  estimates  made. 

| 
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2 . Sample  Problems 
a.  Exact  Bearing  Information. 

Own  ship  tracks  on  course  000°  at  15  knots  for  10  minutes  and 
turns  to  057°  at  15  knots.  Three  bearing  observations  are 
made  on  the  first  leg,  with  the  fourth  observation  made  on 
the  second  leg.  Exact  bearing  information  is  assumed  avail- 
able. The  contact  is  tracked  as  follows: 

Time  Bearing 


1000 

150.0 

1005 

148.2 

1010 

146.9 

1020 

154.7 

Estimate  the  target's  range,  course  and  speed  at  time  1020. 
(Answers:  12,920  yds,  045°,  and  10.0  knots.) 

The  following  additional  observations  are  made: 


Time 

Bearing 

1025 

158.8 

1030 

163.0 

10  35 

167.2 

Estimate  the  target's  range,  course  and  speed  at  time  1035. 
(Answers:  12540  yds,  045°,  10.0  knots.) 
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iSBi 


Initial ize 


000.05  ENT* 

Own  ship's  course  and  speed 

15.(30  iS£<. 

te.ee  eh:: 

Observation  #1  time  and  bearing 

15e.ee  iSt* 

id. 05  eh:: 

Observation  #2  time  and  bearing 

148.2c  83E2 

CARD  I 

card  10.  J0  £„~i 

Observation  #3  time  and  bearing 

14o.  50  bi£.' 

057.00  c.'iT' 

New  own  ship  course  and  speed 

15.00  OSS; 

10.20  ENT.* 

Observation  #4  time  and  bearing 

154.70  6SEC 

iSBE 

Start  estimation 

CSEh 

Continue  estimation 

10.20  o* 

Time  of  estimate 

12501.01  **» 

Estimated  target  range  in  yards 

CARD  2 6SE2 

48.37  *** 

Estimated  target  course  in  degrees 

fiSBE 

5.00  **» 

Estimated  target  speed  in  knots 

10.25  £N7: 

Observation  #5  time  and  bearing 

155.00  oss; 

CARD  1 10. 30  EH i 1 

Observation  #6  time  and  bearing 

t63.ee  tiS£2 

10.35  EHT' 

Observation  #7  time  and  bearing 

107.20  ifpr 

OS  BE 

Start  estimation 

bSE.n 

Continue  estimation 

10.35  *** 

Time  of  revised  estimate 

12540.45  *♦* 

Revised  estimated  target  range 

CARD  2 CS£r 

40.22  *** 

Revised  estimated  target  course 

GSBE 
5.55  *** 


Revised  estimated  target  speed 


Inaccurate  Bearing  Information 

Own  ship  tracks  on  course  000°  at  15  knots  for  10  minutes 
and  turns  to  056°  at  15  knots.  Three  bearing  observations 
are  made  on  the  first  leg,  with  the  fourth  observation  made 
on  the  second  leg.  The  true  target  bearings  have  been 
altered  by  normal  random  variable  with  mean  zero  and  variance 
0.5  degrees  squared.  The  contact  is  tracked  as  follows: 


Time 

Bearing 

1000 

150.0 

1005 

148.7 

1010 

14  6.0 

1020 

155.2 

Estimate  the  target's  range,  course  and  speed  at  time  1020. 
(Answers:  12,920  yds,  045°,  and  10.0  knots.) 

The  following  additional  observations  are  made: 

Bearing 
159.7 
163.0 
166.4 

Estimate  the  target's  range,  course  and  speed  at  time  1035. 
(Answers:  12,740  yds,  045°,  and  10.0  knots.) 


Time 

1025 

1030 

1035 
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i»j£« 

initialize 

8C6.80  El<J‘ 

Own  ship's  course  and  speed 

15.80  «-oE< i 

10.00  ENP 

Observation  #1  time  and  bearing 

250.60  USE 

10.05  £wrr 

Observation  #2  time  and  bearing 

146.20  tS£C 

CARD  1 

10.10  EKTf 

Observation  #3  time  and  bearing 

146.O0  tStC 

056.00  ENT r 

New  own  ship  course  and  speed 

15.00  £o£u 

10. £0  ent; 

Observation  #4  time  and  bearing 

255. £6  bSBJ 

iSBE 

Start  estimation 

6SBn 

Continue  estimation 

10.20  4»* 

Time  of  estimate 

525.".  14  444 

Estimated  target  range  in  yards 

CARD  2 

330.26  444 

Estimated  target  course  in  degrees 

USEE 

83.67  444 

Estimated  target  speed  in  knots 

10.25  EWJ7 

Observation  #5  time  and  bearinq 

15?.  "0  fiSb'2 

10.30  ENT f 

Observation  #6  time  and  bearing 

163.00  6SE.. 

CARD  1 

10.35  EN  T f 

Observation  47  time  and  bearing 

166.40  6662 

6SBE 

1 

Start  estimation 

CSBn 

Continue  estimation 

10.35  444 

Time  of  revised  estimate 

22?45.?0  44* 

Revised  estimated  target  range 

CARD  2 iSEC 

66.64  444 

Revised  estimated  tarqet  course 

• 

08BE 
9.87  **> 


Revised  estimated  tarqet  speed 


3.  Program  Storage  Allocation  and  Listing 
Registers 


RO: 

fcl 

SO  : 

PS1;  PS3; 

PS10 

R1 : 

EW.  - EZ.  tan  B. 

li  l 

SI : 

EAtA  tan2 

B. 

l 

R2  : 

-EW.  tan  B.  + EZ . tan2 
ill 

B . 
l 

S2  : 

EAt?  tan2 

l 

B . 

1 

R3 : 

EW.At.  - EZ.At.  tan  B. 

li  li  i 

S3: 

-EAt2  tan 

l 

B. 

1 

R4  : 

-EW.At.  tan  B.  + EZ.At. 
li  i li 

tan2  B. 

l 

S4  : 

EAt. 

l 

R5: 

ti-l 

S5: 

EAt2 

l 

R6  : 

t . 
l 

S6 : 

E tan  B. 

l 

R7 : 

W. 

1 

SI: 

E tan2  B. 

l 

R8 : 

z . 

1 

S8 : 

EAt^  tan  B 

i 

R9 : 

PS1 1 ; PS12;  v;  At. 

S9  : 

i = N 

RA: 

W 

RB: 

i 

RC: 

A t^ ; PS4;  u 

RD: 

tan  B^ ; PS 7 ; 

A 

V 

RE: 

Z.  tan  B.;  PS13;  X. 

l l 1 

RI : 

i;  PS2 ; PS5 ; 

PS 6 ; PS 8 ; 

All  summations  are  over  the  range  of  i = 1 
PS  denotes  ' partial-sum’ . 
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NOTE 


f • • • t 


Initial  Flag  Status  and  Usp- 
0 : OFF , Unused 
1 : OFF , Used 


2:  OFF,  Unused 
3 : OFF , Used 


User  Control  Keys;  Card  1: 
A: 

B: 

C:  t.  t B.  -»■ 

l l 

D: 

Es  Start  ->■ 


a : 
b: 


Cn  f S, 


c: 

d: 

e:  Initialize 


User  Control  Keys;  Card  2 : 


OCJ 

»L6»: 

L.  4 C * 

002 

sroe 

35  Qi 

002 

ST05 

35  65 

004 

RTH 

24 

005 

tiBLe 

n « 

16  15 

00*' 

CLRG 

16-53 

007 

FZS 

16-51 

006 

CLRG 

16-52 

009 

CLA 

-51 

010 

SF2 

16 

21  02 

011 

RTH 

24 

012 

21 

16  11 

015 

44 

014 

S70B 

35  12 

015 

Ri 

-31 

01c 

STQA 

35  11 

017 

RTH 

24 

016 

tLBLC 

21  12 

019 

ThH 

43 

020 

ST  OB 

35  14 

021 

x:\‘ 

-41 

022 

HNS * 

16  36 

02  3 

ST06 

35  66 

024 

F2° 

16 

23  02 

025 

QSB1 

23  01 

026 

RCL0 

56  60 

027 

- 

-45 

026 

STOC 

«$5  2 J 

029 

2+ 

56 

030 

STOI 

35  46 

031 

Ri 

-31 

032 

RCLC 

56  12 

033 

X* 

53 

034 

X 

-35 

035 

FZS 

16-51 

036 

ST-3 

35-45  02 

037 

RCLD 

36  14 

038 

X 

-25 

039 

ST*2 

35- 

-55  02 

040 

RCLC 

36  13 

041 

RCLD 

36  14 

042 

X* 

53 

043 

X 

-35 

044 

ST+1 

35- 

-55  01 

045 

FZS 

16-51 

046 

RCL6 

36  06 

047 

RCL5 

36  05 

048 

- 

-45 

049 

ENTt 

-21 

050 

EHTt 

-21 

Stores  initial  time,  t^ , for  further  use 

Initialization 

Calculates  and  stores  own  ship  velocity 
components 

Input  observation  and  bearing 
Fills  the  matrix  values 

Branches  to  store  initial  time 


Card  1 


05: 

FCLh 

36  li 

052 

-35 

053 

ST*7 

35-55  6.' 

054 

-31 

055 

RCLG 

36  12 

05* 

A 

-35 

05? 

S7V8 

35-55  OS 

050 

RCL6 

36  06 

059 

ST05 

35  05 

060 

RCL? 

36  07 

061 

S7+1 

35-55  01 

062 

RCLD 

36  14 

067 

X 

-35 

064 

ST-2 

35-45  02 

065 

RCLC 

36  13 

066 

X 

-35 

067 

ST-4 

35-45  04 

068 

RCL  7 

36  07 

069 

RCLC 

36  13 

070 

X 

-35 

071 

sr+3 

35-55  03 

072 

RCL8 

36  0i' 

073 

RCLD 

36  14 

074 

A 

-35 

075 

STOE 

35  15 

076 

ST-1 

35-45  01 

077 

RCLD 

36  14 

078 

X 

-35 

079 

ST+2 

35-55  02 

080 

RCLC 

36  13 

081 

X 

-35 

082 

ST+4 

35-55  04 

083 

RCLE 

36  15 

084 

RCLC 

36  13 

085 

X 

-35 

086 

ST-3 

35-45  02 

087 

RCLI 

36  46 

088 

RTN 

24 

089 

tLBLE 

21  15 

Start  estimate  calculations 

090 

f:s 

16-51 

. 

091 

1 

61 

Changes  the  sign  of  two  matrix  elements 

092 

CHS 

-LL 

093 

ST«6 

35-35  86 

094 

STx8 

35-35  0S 

095 

RCL8 

36  08 

096 

RCL7 

36  07 

(Partial  sums  stored  throughout;  PS) 

097 

> 

-35 

098 

RCL6 

36  06 

099 

RCLI 

36  81 

100 

X 

-35 

84 


— — - 


if: 

PClS 

36  0$ 

15: 

RCLD 

36  1 4 

153 

X 

-35 

154 

- 

-45 

155 

RCL4 

36  o4 

15<? 

KlC 

36  13 

15' 

X 

-35 

15b* 

- 

-45 

15? 

RCLe 

36  06 

It'S 

RCLI 

36  4 c 

161 

x 

-35 

162 

- 

-45 

163 

RCL? 

3b  0$ 

A 

164 

T 

~i.4 

X1 

165 

STQl 

35  15 

166 

1 

01 

Change  sign  of  two  matrix  elements 

167 

CHS 

-22 

16$ 

$Tx« 

35 

-35  c : 

169 

STxS 

35- 

-35  8c 

ire 

prs 

16-51 

171 

RCL6 

36  86 

172 

■>HMS 

16  35 

t . 

173 

FRTX 

-14 

l 

174 

RCL6 

36  06 

175 

RCL8 

36  00 

176 

- 

-4: 

177 

ST  09 

35  89 

17S 

RcLC 

36  13 

179 

X 

-35 

ISO 

RCLE 

36  15 

181 

+ 

.rr 

182 

RCL7 

36  07 

183 

- 

-45 

184 

ys 

53 

18  5 

RCL9 

36  09 

186 

RCLO 

36  14 

187 

V 

-35 

18S 

PC  L! 

36  46 

189 

♦ 

-55 

198 

RCLS 

36  08 

• 

191 

- 

-45 

192 

X* 

53 

193 

♦ 

-55 

194 

195 

;x 

2 

54 

82 

rt  = ( (Xx  + u,  t.  - W.)2  + (Yj  + v,  t.  - Z.)V/2 

19b* 

EEX 

-23 

(in  nautical  miles) 

197 

63 

A 

198 

A 

-35 

R^,  in  yards 

199 

FRTX 

-14 

200 

SF1 

lb 

21  01 

201 

RTt, 

202 

*LBlC 

21  13 

203 

PCLl 

36  13 

204 

FOLD 

36  14 

20  5 

+F 

34 

206 

x:r 

"4  I 

20 7 

x<&> 

16-45 

200 

i£P2 

23  62 

20 9 

FRTX 

-14 

210 

CFl 

16 

22  61 

211 

Rrh 

24 

212 

*LSL2 

21  02 

213 

3 

63 

214 

6 

06 

215 

6 

06 

216 

♦ 

-55 

217 

RTN 

24 

216 

$LBLE 

21  If 

219 

FI? 

16 

23  61 

229 

9T09 

22  69 

221 

R4 

-31 

222 

PRTX 

-14 

223 

SBC 

16-11 

224 

RTN 

•->  .« 

£•»  j 

AD-A068  633 
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E . Mathematical  Analysis 


a.  Assumptions 

All  calculations  use  a rectangular  coordinate  system  as 
defined  below.  Additionally  it  is  assumed  the  following  quantities 
are  accurately  known  (although  the  accuracy  of  the  observed 
bearing  varies) : 

(1)  time 

(2)  target  bearings 

( 3)  observer  course  and  speed 

(4)  observer  initial  position. 

b.  Symbology 

(1)  Observer 

(a)  W:  East-West  position 

(b)  Z:  North-South  position 

(c)  (W^,Z^):  position  at  ith  observation. 

(2)  Target 

(a)  X:  East-West  position 

(b)  Y:  North-South  position 

(c)  (X^,Y^) : position  at  ith  observation 

(d)  u:  East-West  velocity  component 

(e)  v:  North-South  velocity  component. 

(3)  Other 

(a)  : measured  bearing  from  observer  to  target  at  the 

ith  observation. 

(b)  t^:  time  of  the  ith  observation 

(c)  At^:  elapsed  time  between  ith  and  first  observation. 

(d)  i = 1 is  the  initial  observation. 
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c.  Development 

The  geometry  for  a two  leg  TMA  is  shown  in  Figure  1. 

The  target  motion/position  at  any  time  can  be  described  in  terms 
of  its  initial  position  (X1#Y^)f  and  its  velocity  components 
(u  and  v)  , which  are  unknown  and  the  elapsed  time,  At^,  which 
is  known  by:  = X^  + uAt^  and  Y^  = Y^  + vAt^.  Knowledge 

of  the  target  bearing  leads  to  the  following: 

n_i  i _ 1 l l 

tan  B - y"  -~z : - r,  + vAv."--r: 

111  11 

or 

X,  - Y.  tan  B,  + uAt.  - vAt.  tan  B.  = W.  - Z.  tan  B.  (1) 

11  l l l ill  l 


Equation 

(1) 

has  four  known 

variables 

"V  V 

At.  , B 

l 

and  the  four  unknown 

target  variables 

(XL , 

V 

u,  v)  . 

Define 

estimates  for  the 

four  unknowns  as 

A 

Xl* 

V 

U, 

A 

and  v 

and 

define  an  error,  e^,  that  represents  the  errors  due  to  the  use 
of  estimates  in  Equation  (1)  at  each  observation: 


e . = X,  - Y.  tan  B . + uAt . - vAt . tan  B . - W.  + Z . tan  B . . 
ill  ill  ill  l 

The  least  squares  estimates  of  the  variables  X^,  Y^,  u 
and  v are  those  values  which  minimize  the  expression 

n 2 

I e . = E (X.  , Y.  , u,  v) 
i=  1 1 

where  n is  the  number  of  observations. 


94 


s 


Observer  at  first  observation 
time  t^  position  (0,0) 


FIGURE  1.  TMA  Geometry 


To  minimize,  take  the  partial  derivatives  of  E with 
respect  to  each  variable  and  set  each  partial  derivative  equal 
to  zero.  The  resulting  set  of  four  linear  equations  in  four 
unknowns  (Figure  2)  can  be  solved  at  each  bearing  observation 
if  i > 4 and  an  observer  course  or  speed  change  has  occurred 
during  the  TMA  period. 

The  following  target  parameters  can  then  be  calculated 
at  each  successive  observation: 


Target  Range  * \/(X^  - UT)  + (Y^  - Z^) 


Target  Speed 


/T2 ^ 

v u + v 


Target  Course  = arctan(u/v)  with  appropriate  logic 

to  select  the  correct  coordinate 
quadrant . 

When  calculating  these  parameters  appropriate  constants  are 
required  to  insure  proper  units. 


F.  Reference . 

1.  P.  W.  Marzluff  and  R.  C.  Pilcher,  "Basic  Calculator  Methods 
of  Bearings-Only  Tarqet  Motives  Analyses  for  a Moving 
Sensor  (U) . Naval  Postgraduate  School  Thesis,  December  1978. 
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II 

•H 


dt^j  ii 
•H 

I 


Ct^J  II 
•H 


et'-i  it 

Ct^J  II 

Ct-o  II 

C2c — > II 

•H 

•H 

1 

•H 

•H 

1 

•H 

•H 

a 

CQ 

•H 

fN 

•H 

CQ 

C 

C 

a 

3 

<TJ 

fM 

P 

P 

e 

B 

3 

•H 

•H 

43 

4J 

P 

P 

<J 

<3 

H 

H 

H 

H 

C«-~J  II 

c^-j  M 

Ct— J II 

Ct-J  II 

■H 

•H 

•H 

•H 

1 

1 

•H 

CQ 

C 

•H 

3 

(0 

c 

•H 

P 

•H 
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VIII.  NAVAL  GUNFIRE  SUPPORT  GRID  SPOT  CONVERSIONS,  TRUE  WIN D, 


AND  TIME  OF  FLIGHT/MAXIMUM  ORDINATE  COMPUTATIONS  FOR 
5 -INCH/5 4 PROJECTILE  by  LT  Keith  P.  Curtis 

A.  Problem  Statement 

The  success  of  Naval  Gunfire  Support  operations  in  the 
Combat  Information  Center  (CIC)  is  a function  of  rapid  informa- 
tion processing  and  relay.  Specifically,  substantial  error  can 
be  introduced  by  inaccurate  grid  spot  conversions  and  to  a lesser 
degree  by  improper  computations  of  true  wind.  Also,  commence- 
ment of  a fire  mission  can  be  delayed  waiting  for  Time  of  Flight 
(TOF)  and/or  Maximum  Ordinate  (Max  Ord)  information. 

The  inherent  error  of  rapid  calculations  can  be  minimized 
by  the  use  of  the  handheld  programmable  calculator.  This  paper 
addresses  the  use  of  the  Hewlett  Packard  HP-67  to  perform  grid- 
spot  conversions;  compute  true  wind,  TOF,  and  Max  Ord. 


B.  Operational  Analysis 

The  objective  of  this  program  is  to  provide  a orie-card 
program  to  accommodate  the  following: 


1.  Correct  for  magnetic  variance  for  a geographic  area. 

2.  Accept  the  Observer  Target  Line  (OTL)  in  either  mils 
magnetic  or  degree  magnetic. 

3.  Perform  precise  grid  spot  conversion. 

4.  Provide  Time  of  Flight  information. 
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5.  Provide  Maximum  Ordinate  information. 


6.  Compute  true  wind. 

C.  Computational  Algorithm 

1.  Enter  magnetic  variance. 

2.  Enter  OTL  either  in  mils  magnetic  or  degree  magnetic. 

3.  Enter  observer  "spots":  left-right,  add-drop  (in  yards). 

4.  Convert  spots  to  East-West,  North-South. 

5.  Enter  range  of  shot  and  compute  TOF  and  max  ord. 

6.  Enter  own  ship  course  and  speed,  and  relative  wind  to 
compute  true  wind. 

7.  Repeat  Steps  3 through  6 as  necessary. 
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A. 


HP-67  Calculator  Proqram 


. User  Instructions 


Step 

Instruction 

Input 

Key(s) 

1. 

Read  program  card  (both  sides) 

2. 

Enter  magnetic  variation  (+  for  East, 
- for  West) . If  no  variation  is 
entered,  0 is  used. 

mag. var. 

f b 

3. 

Enter  Observer  Target  Line 
in  either: 

a. 

degrees  magnetic,  or 

A 

b. 

mils  magnetic 

f a 

4 a. 

Enter  Left/Right  spot 
(-  for  left,  + for  right) 

L/R 

t 

b. 

Enter  Add/ Drop  spot 
(-  for  drop,  + for  add) 

A/D 

B 

c. 

Display  E/W  spot  (+  E,  - W) 

d. 

Display  N/S  spot  (+  N,  - S) 

Optional:  Recover  the  E/W  spot 

h x X y 

e. 

Repeat  Step  4 as  required 
until  OTL  changes 

5. 

Compute  TOF  and  max  ord 

a. 

Enter  target  range,  display: 

range 

C 

b. 

Time  of  Flight 

c. 

Max  Ord 

Optimal:  Recover  TOF 

h R t 

6. 

Compute  true  wind. 

a. 

Enter  own  ship' s course  and 
speed  in  the  form  SS.CCC 
where  SS  is  the  speed  in 
integer  knots  and  CCC  is  a 
three-digit  course 

SS.CCC 

+ 

b. 

Enter  relative  wind 

SS.CCC 

D 

c. 

True  wind  is  displayed 

E/W  spot 
N/S  spot 


TOF 

Max  ord 


SS.CCC 


Note:  Use  of  the  keys  [c]  or  [d]  do  not  require  an  OTL  input. 
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Sample  Problem 

Own  ship: 

0 35T 

12KTS 

Relative  wind: 

225R 

6KTS 

Magnetic  variance: 

7E 

OTL : 

1930 

mils  mag 

Spot : 

Left 

250,  Drop  150 

Range : 

9600 

yds 

West  27, 

North  290 

TOF: 

16  seconds 

Max  Ord: 

1050ft 

True  Wind: 

2 30T  17KTS 

7. 

6661 

193 d. 

•S£;. 

116. 

*** 

-250. 

ENT* 

-150. 

-27. 

»«« 

290. 

*»* 

9600. 

csg: 

16. 

m 

1056. 

«*• 

12.035 

ENT* 

6.225 

css: 

17.230 

♦ »* 

3.  Program 

Storage  Allocation 

and  Listings 

Registers 

R0; 

OTL  °T 

SO: 

RA: 

log1Qrange 

Rl : - 

SI: 

RB: 

Ship's  heading 

R2 : 

S2  : 

RC: 

L-R  spot 

R3 : 

S3: 

RD: 

A-D  spot 

R4  : 

TOF/ 

S4  : 

IX 

RE: 

“ Max  Ord 

R5: 

Coef  f 

S5  : 

RI : 

Control 

R6  : 

S6  : 

ZY 

R7 : 

S7: 

R8 : - 

S8: 

R9: 

mag. var . 

S9 : 

Initial  Flag  Status  and  Use 
0:  OFF,  Unused 

1:  OFF,  Unused 


2:  OFF,  Unused 

3:  OFF,  Unused 


User 

■ Control  Keys 

A: 

OTL  (degrees) 

a:  OTL  (mils) 

B: 

L-R  spot,  A-D  spot 

b:  mag. var. 

C: 

Range 

c : 

D: 

O/S  c/S  rel  wind 

d : 

E: 

e : 
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E . Geometric/Mathematical  Analysis 
1 • Grid-Spot  Conversion 

The  conversion  of  grid  spots  oriented  to  an  Observer 
Target  Line  to  an  East-West,  North-South  orientation  for  input 
to  a shipboard  GFCS  can  be  accomplished  by  a rotation  of  the 
OTL  counterclockwise  to  000  degree  True  after  the  OTL  has  been 
corrected  to  true  bearing  (Figure  1) . 

The  rotation  formulas  are 


or 


cos  0 sin  0 

X 

X ' 

-sin  0 cos  0_ 

. y . 

_ y'  . 

x'  x cos  6 + y sin  0 and  y'  = -x  sin  0 + y cos  0, 


where 


6 = OTL  °T 
x = L/R  Spot, 
y = A/D  Spot, 
x’  = E-W  Spot, 
y'  = N-S  Spot 
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2 . TOF  and  Maximum  Ordinate 

This  information  is  tabulated  for  the  5"/54  Projectile 
in  BuOrd  Publication  OP1182  (Range  Table  for  5"/54) . To  find 
an  equation  which  would  best  fit  the  tabulated  data,  a log 
transformation  was  made  followed  by  a parabolic  curve  fit. 

The  time  of  Flight  Equation  is 


f(x)  = 10 


(0.1083  + 0.8505  log  x + 0.2831  (log  x)  ) 


This  equation  generates  solutions  within  1 second  of  the  tabulated 
values  for  ranges  under  22,000  yds.  The  same  approach  was  made 
in  determining  an  equation  for  maximum  ordinate. 

The  Maximum  Ordinate  Equation  is 


f (x)  = 10 


(0.  7775  + 2 . 16661og x -0.2071 (loq  x)  2 +0.261 7 (log  x)  3 + 0.0721 (log  x) 4) 


This  equation  generates  solutions  within  55ft  of  the  tabulated 
values  for  ranges  less  than  19,000  yds.  (Within  10ft  for  range 
less  than  13,000  yds,  and  with  1 ft  for  ranges  less  than  7000  yds.) 

The  algorithm  uses  a nested  polynomial  to  preserve 
accuracy  in  the  calculation  of  the  exponent.  For  example. 

Time  of  Flight  is  computed  as  follows: 


f(x)  = 10<<0'2831  log  x + 0*  8505) log  x + 0.1083) 
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3.  True  Wind. 


This  computation  is  simple  vector  arithmetic.  Relative 
wind  is  converted  to  Apparent  wind  by  adding  the  ship's  heading 
to  the  relative  bearing,  then  converted  to  rectangular  coordi- 
nates. Own  ship's  course  and  speed  are  then  converted  to 
rectangular  coordinates  and  subtracted  from  apparent  wind  vector. 
The  result  is  true  wind  which  is  converted  to  polar  coordinates 
(Figure  2) . 


N 


0 = relative  wind  bearing 
Y = apparent  wind  bearing 

FIGURE  2.  Wind  vectors 
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1 

Step 

Instruction 

Input 

Key (s) 

Output 

1. 

Read  program  card 

2. 

Initialize 

f e 

1.00 

3.  a. 

Sound  velocity  (m/sec) 

c0 

t 

b. 

W»ter  depth  (m) 

d 

t 

c. 

Source  frequency  (Hz) 

f 

4. 

Bottom  type:  Esther 

a. 

Rigid  Bottom,  or 

none 

A 

b. 

Pressure  release  bottom 

B 

5. 

Output  sequence 

a. 

Mode  number 

n 

b. 

Cutoff  frequency  for  mode  n 

fC 

c. 

Group  velocity  for  mode  n 

cg 

d. 

Phase  velocity  for  mode  n 

°P 

e . 

Display  mode  number 

n 

6. 

Continue  from  Step  4a  or  4b 
depending  upon  original 
bottom  type . 

1 

This  process  continues  until 
the  highest  mode  that  will 
propagate  for  the  given  con- 
ditions has  been  displayed 
when  Q or  [b]  is  pressed 
and  the  next  mode  will  not 
propagate  then  the  program 
displays  the  mode  number  of 
the  last  mode  that  will 
propagate. 

2 . Sample  Problems 
Example  1 . 

Cq  «*  1500  m/sec 
d = 15  m 
f = 150  Hz 


Rigid  bottom 


Initializes  Program 

1500.00  Ei,:: 

15 .00  EH7‘ 

1 5LkC0  OCZn 

1.60  **» 

Flashes  mode  number  (n) 

25.66  **♦ 

Flashes  fc  for  n * 1 

1479.62  »** 

Flashes  Cg  for  n * 1 

1521.20  »»♦ 

Flashes  Cp  for  n = 1 (Stops,  displaying  mode  no.) 

6SBA 

2.66  »*< 

Flashes  n 

75.06  •** 

Flashes  fc  for  n * 2 

1299. 64  *** 

Flashes  Cg  for  n = 2 

1732.65  *»* 

Flashes  C..  for  n - 2 (Stops,  displaying  mode  no.) 

CS86 

3.06  *»* 

Flashes  n 

125.66  »«4 

Flashes  fr  for  n * 3 

829.16  **» 

Flashes  Cg  for  n = 3 

2713.66  *»* 

6S6(i 

Flashes  Cp  for  n - 3 (Stops,  displaying  mode  no.) 

3.66  *** 

Mode  4 will  not  propagate  under 

CSB  h 

these  conditions  so  regardless  of  how  many  times  [a] 

3.66  **» 

is  pressed,  mode  3 is  displayed. 

Example  2 


Cq  « 1500  m/sec  Pressure  Release  Bottom 

d » 15  m 
f = 150  Hz 


CSEe 

Initializes  Program 

1500.00  en:: 

15.00  ENT  I 

156.00  CSB3 

1.00  »«* 

Flashes  mode  no. 

50.00  **» 

Flashes  for  mode  1 

1414.31  »»• 

Flashes  Cg  for  mode  1 

1590.95  •** 

Flashes  CD  for  mode  1 (Stops,  displaying  mode  no.) 

CSBB  ' 

2.00  **» 

Flashes  mode  no. 

100.00  •»» 

Flashes  fc  for  mode  2 

1110.05  *»♦ 

Flashes  Cg  for  mode  2 

2012.40  M# 
CSBF 

Flashes  Cp  for  mode  2 (Stops,  displaying  mode  no.) 

2.00  *♦* 

Only  two  modes  will  propagate 

2.00  »»* 

NOTE:  If  such  conditions  exist  such  that  no  modes  will  propagate, 

0.00  is  displayed. 


The  program  stops  between  modes  and  requires  the  user  to  initiate 
the  next  mode  in  order  to  allow  the  user  sufficient  time  to  write 


down  the  information  presented. 


3 . Program  Storage  Allocations  and  Program  Listim 
Registers : 


R0;  f 

SO: 

RA: 

Rl:  d 

SI : 

RB: 

R2:  CQ 

S2: 

RC: 

R3 : n 

S3: 

RD: 

R4:  fc 

S4  : 

RE: 

R5:  /l  - <fc/f)2 

S5  : 

RI  : 

R6: 

S6  : 

R7 : 

S7: 

R8 : 

S8: 

R9  : 

S9: 

Initial  Flag  Status  and  Uses 
0:  OFF,  Unused 

Is  OFF,  Unused 


2.  OFF,  Unused 
3:  OFF,  Unused 


User  Control  Keys 
A:  Rigid  Bottom 

B:  Pressure  Release  Bottom 

C: 

D: 

E: 


. 


e«: 

• LSI.* 

21 

10  15 

Initializes  Program 

00’ 

1 

01 

003 

ST03 

35  03 

n « 1 

004 

SF2 

16 

21  02 

Controls  storage  of  inputs  on  first  pass 

00? 

It'S 

51 

Input  parameters  at  this  stop 

00C 

tiBLO 

21  11 

Case  I — Rigid  Bottom 

007 

F2~ 

16 

23  02 

90S 

C$61 

23  01 

Stores  inputs  on  first  pass 

009 

RCL3 

36  02 

010 

• 

-62 

Oil 

5 

65 

01 2 

X 

-35 

013 

• 

-62 

014 

£ 

02 

Calculates  f 

015 

5 

05 

C 

016 

- 

-45 

C-(2n-l) 

017 

RCL2 

3c  02 

fC  “ 4d 

01S 

X 

-35 

019 

RCL1 

36  01 

020 

* 

-24 

j 

021 

r:l0 

36  00 

022 

x?y 

-41 

027 

x*v? 

16-32 

Is  f > f 

024 

x>v? 

16-34 

C — 

025 

CJ04 

22  04 

Yes — display  last  mode  # and  stop 

020 

ST04 

35  04 

No — continue 

027 

RCL3 

3c  03 

022 

FR  TX 

-14 

Flash  new  mode  # 

025 

RCL4 

36  04 

030 

FR7X 

-14 

Flash  fc 

071 

CS  62 

23  02 

Computes  Cg  and  Cp 

032 

*LBL3 

21  03 

Increments  mode  no.  for  next  iteration 

033 

1 

01 

034 

ST*3 

35 

-55  03 

035 

RCL3 

3c  03 

03c 

1 

01 

037 

- 

-45 

03S 

R S 

51 

.Stops,  displays  mode  ft  of  tun  just  completed 

039 

$lblb 

21  12 

Case  II--Pressure  Release  Bottom 

040 

F2? 

16 

23  62 

041 

CSfl 

23  01 

-Stores  inputs  on  first  pass 

042 

RCL3 

36  63 

047 

RCL2 

36  02 

044 

X 

-35 

Calculates  f 

045 

RCL1 

36  61 

C 

046 

A 

-24 

* 

047 

2 

02 

C (2n) 

040 

A 

-24 

■c 

) 

■ 

045 

RCLO 

36  06 

050 

x;v 

-41 

051 

xpv* 

16-32 

052 

X>Y? 

16-34 

053 

6104 

22  04 

054 

ST04 

35  04 

055 

RCL3 

36  03 

056 

FR1X 

-14 

05? 

FCL  4 

36  6-i 

058 

FF1X 

-14 

05? 

6SB2 

23  02 

06t» 

61 83 

22  03 

061 

»L6ul 

21  01 

063 

ST80 

35  66 

063 

R4 

-31 

064 

ST01 

35  01 

065 

R4 

-31 

066 

ST  62 

35  C-2 

06? 

RTH 

24 

068 

*LBL2 

21  02 

06? 

RCL4 

36  04 

0?O 

RC  L0 

36  0O 

071 

? 

072 

X* 

53 

073 

CHS 

074 

1 

01 

075 

♦ 

-55 

076 

;x 

54 

077 

ST05 

35  05 

076 

RCL2 

36  02 

07? 

X 

-35 

080 

PRTX 

-14 

081 

RCL2 

36  62 

082 

RCL5 

36  65 

083 

* 

_ 

084 

PR  TX 

-14 

085 

RTH 

24 

086 

»LBL4 

21  04 

08' 

RCL3 

36  07 

088 

1 

01 

08? 

- 

-45 

0?0 

R$ 

51 

Is  fc  > f 

Display  last  mode  no. 

Yes — and  stop 
No — continue 

Flash  new  mode  no. 

Flash  fc 

Computes  Cg  and  Cp  increments  mode  # for  next 
iteration 

Subroutine  1 — Stores  inputs  on  first  pass  only 


Subrounte  2 — Calculates  C and  C 

g P 


C 

g 


Flash  Cg 

CP  - c„A/i  - 

Tlash  C 

P 


Allows  the  last  mode  number  to  be  displayed 
after  each  iteration 
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E.  Mathematical  Analysis 

Consider  a shallow  water  ocean  model  with  depth  d. 
According  to  Normal  Mode  Theory  the  sound  pressure  at  any  point 
can  be  determined  by  the  solution  of  the  wave  equation  and  is 
given  by 


pn  - An  sin  Knz  z e xp [ i ( u> t - K^x)  J 

Sound  will  propagate  at  different  angles  and  this  is 
what  gives  rise  to  the  various  modes.  As  0,  in  Figure  1 
approaches  zero,  the  sound  will  not  propagate  because  it  is  merely 
bouncing  up  and  down  off  the  surface  and  bottom  (no  x-direction 
of  travel) . This  determines  the  cutoff  frequency  for  mode  n. 

The  cutoff  frequency  is  dependent  upon  the  boundary 
conditions  at  the  surface  and  bottom  because  the  solution  of 
the  wave  equation  is  dependent  upon  the  boundary  conditions. 

The  surface  is  always  considered  to  be  a pressure  release 
boundary.  As  the  mode  number  increases  the  cutoff  frequency 
for  that  mode  is  higher  also.  When  the  cutoff  frequency  exceeds 
the  frequency  of  the  source  then  that  mode  will  not  propagate, 
nor  will  higher  modes. 
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Case  I:  Rigid  Bottom 


For  the  rigid  bottom,  the  boundary  condition  yields  a 


(2n-l) * 


n * 1,2,3  (mode  #) 


c _ 0 ( 2n-l ) n 

fC  ~ Y~  ~ 73 * 


Case  II:  Pressure  Release  Bottom 


For  the  pressure  release  bottom,  the  boundary  condition 


yields 


n = 1,2,3,...  (mode  #) 


17  IT 


The  remaining  values  calculated  are  the  group  velocity 


Cg  = CQ  cos 


and  the  phase  velocity  C = C./cos  $ , where 

p u 


2 1/2 

cos  $ = [1  - (fc/f)  ] / , fc  is  the  cutoff  frequency,  f is 
the  source  frequency,  and  Cq  is  the  velocity  of  sound.  The 
geometry  of  these  quantities  is  shown  in  Figure  2. 
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NORMAL  MODE  TRANSMISSION  LOSSES  by  LT  Michael  D.  Clary 


A.  Problem  Statement 

As  an  alternative  to  Ray  Theory,  the  use  of  Normal  Mode 
Theory  provides  a more  exact  approach  to  the  solution  of  trans- 
mission loss  problems.  Given  source  frequency  and  depth,  range 
to  and  depth  of  receiver,  and  modal  values  for  the  speed  of 
sound  and  the  absorption  coefficient,  the  program  user  may 
either  specify  the  effective  pressure  amplitude  of  the  source 
at  one  meter  and  solve  for  modal  effective  pressure  amplitudes 
and  the  resulting  transmission  losses,  or  if  the  modal  pressures 
are  input,  the  value  of  the  one  meter  effective  source  pressure 
amplitude  can  be  obtained. 

B . Operational  Analysis 

Since  this  program  is  used  primarily  as  a theoretical 
problem  solver  based  on  give  data,  no  operational  analysis  is 
provided. 

C.  Computational  Algorithm  HP-67 

1.  Input  source  frequency,  depth,  receiver  depth,  effective 
pressure  amplitude  at  1 meter  (source)  or  effective  pressure 
amplitude  of  given  model. 

2.  Input  absorption  coefficient  for  given  mode,  range  between 
source  and  receiver,  and  speed  of  sound  for  the  mode. 

3.  Compute  sums  required  for  transmission  loss  equation. 

4.  Compute  and  output  coherent  and/or  incoherent  transmission 
losses . 
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HP-67  Calculator  Program 
1 . User  Instructions 


Step 

Instruction 

Input 

Keys (s) 

Output 

1. 

Load  side  one  and  two  of 
program  card 

2. 

Program  the  eigenfunction 

for  z (z) , z (z„)  at  LBL3 
n n 0 

in  the  W/PGRM  MODE. 

(Assume  value  of  z and 

z o will  be  in  the  x-register 
upon  initialization  of  each 

computation.)  The  last  two 
program  steps  must  be  STO(i) , 
h RTN.  Return  calculator  to 

RUN  MODE 

3. 

Enter  conputational  values 

a. 

frequency  of  source  in  Hz 

f 

ENTER 

f 

b. 

depth  of  source  in  meters 

z0 

ENTER 

z0 

c. 

depth  of  receiver  in  meters 

Z 

ENTER 

z 

d. 

effective  pressure  amplitude 
at  1 m in  pPa 

P(l) 

A 

f 

or  effective  pressure  amplitude 
of  given  value 

-lpJ 

A 

f 

enter  as  -1 P 1 

1 n’ 

e . 

absorption  coefficient  for 

an 

ENTER 

a 

the  mode 

n 

n 

f. 

range  from  source  in  meters 

r 

ENTER 

g- 

speed  of  sound  for  the 
mode  in  meters/sec 

C 

n 

B 

|P  I or  P<  1) 

1 n1 

Output  will  be  either  the 
mode  value  for  the  effective 
pressure  amplitude  or  the 
value  of  the  effective  pres- 
sure amplitude  at  1 m.  All 
values  are  in  ypa. 
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Step 

Instruction 

Input 

Key (s) 

Output 

4. 

Sum  the  values  of 

C 

n 

P sin(27rf  )/C 
n r n 

. 

and 

■ 

P cos  ( 2v  f ) /C 
n r n 

5. 

2 

Sum  the  value  of  P 

n 

D 

p2 

n 

6. 

Return  to  Step  2 for 
additional  mode  computations. 
Information  for  Steps  3a, 
b,c  and  d need  not  be  re- 
entered. If  all  mode  com- 
putations are  completed, 
go  to  7. 

B 

Compute  transmission  loss 
(in  dB) 

for  coherence 

E 

T.L. 

A 

for  incoherence 

f 

e 

T.L. 

H 

For  new  case  go  to  Step  1 
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2 . Sample  Problem 

For  a nearly  isovelocity  layer  of  water  50  m deep 
overlying  a silt  bottom  rich  in  decaying  organic  matter,  the 
following  are  found  to  be  good  approximate  values  for  a source 
frequency  of  50  Hz: 


Zn(z) 


f (cutoff) 
n 


0.2 

sin  (0 . 050z) 

1550m/s 

25  Hz 

1 

x 10-4/m 

0.2 

sin (0 . 105z) 

1630 

33 

2 

x 10~4 

0.2 

sin (0. 189z) 

2110 

45 

4 

x 10~4 

not 

excited 

— 

55 

~ — — 

Evaluate  the  effective  amplitudes  I I at  a receiver 
at  a range  of  2 km  and  at  a depth  of  20  m.  Determine  both 
transmission  loss  values. 

SOLUTION:  (given  that  the  source  is  at  a depth  of  50m) 

1.  Enter  first  eiqen function  Z (z)  under  label  3. 

n 

Value  of  z (or  z^)  will  be  in  the  X-register. 

2.  Enter  f,  zQ,  z and  P ( 1 ) = 107  MPa  for  n = 1 and 
compute  P^ . 

3.  After  recording  value  P^,  sum  the  values  for  trans- 
mission loss  calculations. 

4.  Repeat  Steps  1-3  for  all  necessary  values  of  n (three 
for  the  given  example) . 

5.  When  all  P 's  have  been  obtained,  calculate  transmission 

n 

loss,  both  coherent  and  incoherent  cases. 

6.  Numerical  answers: 

4 . . 4 


P3  = 0.39  x 10  MPa 

TL(coherence)  = 53  dB 
TL( incoherence)  = 50  dB 


MPa 

|P1  1 

= 2.05 

x 104 

MPa 

MPa 

1 P2  1 

= 2.54 

x 104 

MPa 

MPa 

lp3l 

= 0.39 

o 

X 

MPa 

Keystroke  Sequence  for  Sample  Problem  a 


b.  Using  the  same  table  values  and  the  computed  values  for 


each  P , verify  that  the  value  for  P ( 1 > is  in  fact 
n 

10 7 v Pa . 

1.  Enter  first  eigenfunction  Z^(z)  under  label  3. 

2.  Enter  computational  values  for  n = 1 and  compute 
Note  that  the  values  of  P must  be  entered  as  the 
negative  value  of  the  absolute. 

3.  Calculate  P(l).  Repeat  Steps  1-3  for  each  value  of 

4.  Numerical  answers: 


(Pi1 

= 2.05 

4 

X 10* 

ViPa 

nil 

1 P2 1 

= 2.54 

x 104 

Vi  Pa 

Pd) 

1 P3  1 

« 0. 39 

x 103 

Vi  Pa 

p (i ) 

Keystroke  Sequence  for  Sample  Problem  h. 


erf: 

50.00  Eft'* 

50.00  £.S  7 * 

20.00  Eft 7* 

-2.  05*04 

1. -04  FN7.“ 

2.  +03  Eft7 : 
1550.00 

?.  »»* 
ero: 

50.  EST? 
50.  ENT? 
20.  ENT? 
-2. 54*04  Ofto 
2.-04  ENT? 

2.  *03  EN T ? 
1630.  i 7666 

1.00*07  ♦ »* 

tros 

50.  est; 
50.  ENT? 
20.  EST! 
•0.30*03  k'io 
4.-04  ENJl 

3.  *03  ENT? 
2110.  i3£f 

>.94*00  »»* 


= 9.98  x 10*  « 10 7 
= 1 0 u Pa 

= 9.94  x lob  « io7 


P ( U . 


n . 


; .Pa 


viPa 
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3.  Program  Storage  Allocation  and  Program  Listirn 


Reg 

isters : 

RO: 

f 

SO: 

Zpn 

R1 : 

zo 

SI: 

R2: 

z 

S2 : 

R3: 

Cn 

S3: 

R4  : 

a 

n 

S4  : 

T.  sin  (2nfr) /Cn 

R5: 

r 

S5: 

R6  : 

P(l)  or  - 1 Pn 1 

S6: 

E cos  ( 2tt  f r)  /Cn 

R7j 

S7: 

Z cos2  ( 2tt  f r ) /CR 

R8 : 

S8 : 

£ [sin(2nfr)/Cn] 

R9  : 

S9: 

n 

RA: 

P or  P(l) 

n 

RB : 

(27Tfr)/Cn  and 

COS  ( 2 Tlf  ) /C 
r n 

RC: 

sin(2Trfr)  /Cn 

RD: 

Coherent  TL 

RE: 

Incoherent  TL 

RI : 

Scratch 

Initial  Flag  Status  and  Use:  OFF,  Unused 
Trig  Mode:  RAD 


001 

»Ltl0 

cl  11 

002 

K00 

16-22 

oo? 

CLR3 

16-5? 

C84 

p:  s 

16-51 

005 

cue 

16-5? 

m 

STl>6 

35  06 

PfT 

Ri 

-31 

063 

ST  02 

35  02 

009 

Ri 

-31 

oio 

ST  01 

35  01 

ou 

Ri 

-31 

012 

S TOO 

35  00 

013 

RTU 

24 

014 

$LBl  6 

21  12 

615 

ST  03 

35  0? 

016 

Ri 

-31 

01? 

ST05 

35  05 

01  s 

Ri 

-31v 

01 9 

ST04 

35  04 

020 

RCL6 

36  06 

021 

X<0n 

16-45 

022 

CT  01 

22  01 

02? 

1 

01 

024 

srn? 

35  09 

025 

*LFLi 

21  02 

0±6 

7 

0? 

02? 

ST  01 

35  4 c 

02S 

RCL1 

36  01 

029 

bi  t 3 

23  65 

030 

6 

06 

031 

STC1 

35  46 

032 

RCl2 

36  o? 

033 

6SB3 

23  C? 

034 

RCL? 

36  02 

03  5 

X<0° 

16-45 

036 

RTN 

V4 

£.** 

03? 

net  4 

21  04 

039 

RlL? 

36  0? 

039 

RCL  5 

36  65 

040 

fX 

r * 

041 

7 

•£4 

042 

r:ls 

36  06 

043 

t 

-35 

044 

RCL  3 

36  65 

045 

RCLO 

36  C3 

046 

7 

-34 

047 

:x 

54 

04$ 

y 

-35 

049 

KCL6 

36  06 

650 

x 

-35 

051 

RrL  4 

36  04 

853 

RCl5 

36  05 

05? 

» 

-35 

054 

CHS 

055 

«• 

33 

Pr>6 

\ 

-35 

057 

35  11 

Enter 

values  for 

f,  zQ, 

Z,  and  either 

P(l) 

or  -|  P | . 

1 n1 

Clears 

all  registers, 

sets 

radian  mode. 

Enter 

values  for 

a , r, 
n 

and  c 

n 

Ctetermine  calculations  required. 

1 is  stored  in  FT  if  P is 
9 n 

beiny  calculated 


11 5 

w 

-35 

Jit 

2* 

56 

u: 

PTH 

7 « 

«.  •» 

ii8 

tLPLO 

21  14 

u? 

RCLO 

36  11 

120 

X* 

53 

121 

p;s 

16-51 

122 

ST+6 

35-55  02 

123 

f:s 

16-51 

124 

RTH 

<;4 

125 

fLBLE 

21  15 

126 

p:  s 

16-51 

! 127 

RCL4 

36  64 

\ 128 

X* 

55 

129 

RCL6 

36  06 

130 

X* 

52 

1 31 

+ 

.cr 

132 

,rX 

54 

PUS 

16-51 

i 

134 

RCL6 

36  06 

135 

r 

*■>  4 

136 

LOb 

16  35 

137 

2 

62 

138 

0 

66 

13? 

X 

-35 

140 

CHS 

~c2 

141 

FIX 

-11 

142 

D8F0 

-63  66 

143 

STOO 

35  14 

144 

RTH 

1 4 

145 

* LBLe 

21  16  15 

146 

f:s 

16-51 

147 

RCLO 

36  00 

148 

f:s 

16-51 

148 

fX 

54 

158 

Rll.6 

36  66 

151 

T 

-24 

152 

LOG 

16  32 

153 

02 

154 

0 

00 

155 

X 

-35 

158 

CHS 

157 

FIX 

-n  j 

158 

P8P0 

-63  06 

15? 

STOE 

35  15 

160 

RTH 

24 

161 

R 8 

51 

• 

Computes 

p2 

n 

Computes 

T.L. 

assuming 

coherence. 

Computes 

T.L. 

assuming 

incoherence . 

E.  Mathematical  Analysis 

The  following  four  formulas  are  the  basis  for  the  mathe- 
matical computations  (Ref.  1) . 


(C-  Z (z  ) 

1.  P = P(l)  \ -Jr-  — Z (z)  exp(-a  r) 

n v f i-  n n 


2. 


Pil)  = 


P Jr 
n 


/ (Cn/f ) Z^ ( Zq ) Zn(z)  exp(-anr) 


3.  Coherent  transmission  loss 


TL  = -20  log 


2f 

I P„  sin  -g 


n i 


' / 2f 

+ I p cos  — ^ 
\ n C. 


n i 


1/2 


'P(l) 


4. 


Incoherent  transmission  loss 


TL  = -20 log 


(E  P^)1/2/P(l) 


F.  Reference 


1.  A.  B.  Coppens  and  J.  V.  Sanders,  "An  Introduction  to  the 
Sonar  Equations  with  Applications,"  Technical  Report 
NPS-61Sd76071 , July  1976,  Naval  Postgraduate  School, 
Monterey,  CA  93940 
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XI.  GOLDEN  SECTION  SEARCH  by  LT  J.  K.  McDermott 


A.  Problem  Statement 

The  minimum  value  of  a unimodal  function  of  one  variable, 
f(x),  for  a specified  interval  is  determined  by  utilizing  Golden 
Section  Search  Techniques,  i.e. 

Minimize  f(x) 

Subject  to  x € I 

where  I = [a,b]  is  a closed  interval  in  E^  space. 


B . Operational  Analysis 

Golden  Section  Search  is  a specific  type  of  interval 
of  uncertainty  (IOU)  method  of  single  variable  optimization 
which  requires  the  selection  of  a specific  interval.  Once 
the  interval  has  been  selected,  the  program  locates  the  value 
of  x which  will  minimize  a unimodal  function  f(x)  being 
evaluated  within  this  specific  interval.  If  a different 
interval  is  selected,  a different  x with  a correspondingly 
different  minimum  functional  value  may  be  obtained. 

Golden  Section  Search  locates  a local  minimum  and  not 
the  global  minimum.  A function  which  is  not  unimodal  over  the 
specific  interval  may  produce  an  x value  which  does  not  provide 
the  minimum  (global)  functional  value  within  the  interval. 

The  behavior  ior  non-unimodal  functions  is  not  predictable. 
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PRECEDING  PjuSK  ai  tif 


The  functions  which  may  be  evaluated  are  limited  to 
some  degree  by  (i)  the  number  of  program  steps  available  for 
user  supplied  function  program  listing  (139-224),  (ii)  functions 
of  one  variable  preferably  unimodal  over  the  IOU  to  avoid 
ambiguity,  and  (iii)  user's  programming  capability  and 
imagination . 

Golden  Section  Search  was  originally  suggested  by 
J.  Kiefer,  "Sequential  Minimax  Search  for  a Maximum,"  Proc. 

Amer.  Math.  Soc. , 4,  no.  3,  June  1953,  pp.  502-506.  The  name 
traces  back  to  Euclid’s  discovery  that  it  is  possible  to  divide 
any  given  line  segment  into  two  parts  such  that  the  ratio  of 
the  whole  to  the  larger  part  equals  the  ratio  of  the  larger 
part  to  the  smaller.  The  division  of  a line  in  this  manner 
came  to  be  known  as  the  Golden  Section,  both  because  it  has 
several  rather  interesting  geometric  and  numerical  properties 
and  because  the  proportions  of  the  two  parts  seem  pleasing  to 
the  eye. 

The  author  (programmer)  is  indebted  to  Professor  J.  K. 
Hartman  of  the  Naval  Postgraduate  School  whose  lectures  and 
class  notes  form  the  bases  of  this  HP-67/97  calculator  program. 
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C . Computational  Algorithm 

Basic  IOU  Algorithm  Structure  (GSS) 

1.  Given  initial  IOU  I = [a,b]  and  function  f (x) . Let  K 
be  a function  evaluation,  iteration  counter. 

2.  Compute  initial  X^  as 

= a + o (b-a)  with  a = (3  - /5)/2 
Set  Ix  = 1 and  K = 2. 

3.  At  iteration  K,  interval  IR-1  resulting  from  previous 
iteration  contains  best  point  (one  producing  smaller  function 
value)  thus  far  and  its  relative  position  is  o or  l-o. 

Place  new  point  (X„)  symmetrically: 

l\ 

XK  = ENDL  + ENDR  - X^ 

where  Iv  . = [ENDL,  ENDR]  and  X provided  the  smaller 

function  value  between  previous  two  evaluated  points. 

4 . Compute  f (XR) . 

5.  Shorten  IOU  to  IK  <-  with  length  LR  < L^_^  from 

information  ftx^  provides.  Set  K = K+l  and  go  to 
step  3 for  the  next  iteration. 

STOPPING  RULE: 

Stop  when  either  K = NMAX  (present  number  of  function 
iterations)  or  when  L^  £ RIOU  (preset  required  interval  of 
uncertainty  length). 

1 39 

_ i 

BhLia. ~ — - - - — -■  - . — 


HP-67  Computational  Algorithm 

1.  Input  user  supplied  function  program  listing  in  available 

program  steps  139  through  224. 

I 

2.  Input  left  endpoint  of  interval  of  uncertainty  (ENDL) . 

3.  Input  right  endpoint  of  interval  of  uncertainty  (ENDR) . 

4.  Input  required  length  of  the  final  interval  of  uncertainty 
( RIOU) . 

5.  Input  maximum  number  of  function  evaluations  desired  (NMAX) 

6.  Output  final  interval  of  uncertainty  [ENDL, ENDR] . 

7.  Output  minimum  (local)  function  value  in  interval. 

8.  Output  X value  that  produces  minimum  function  value. 


. HP-67/97  Calculator  Program 


User  Instructions 


Instruction 


Enter  program  card 
Select  GTO  f e 

Slide  W/PRGM-RUN  switch  to 
W/PRGM 


Slide  W/PRGM-RUN  switch  to 
RUN 

Enter  left  endpoint  of 
interval  of  uncertainty 

Enter  right  endpoint  of 
interval  of  uncertainty 

Enter  required  length  of 
interval  of  uncertainty 

Enter  maximum  number  of 
function  evaluations  to 
be  performed 

Compute  final  interval  of 
uncertainty,  minimum  function 
value,  and  corresponding  X 
value 

ENDL  displayed  when  compu- 
tation complete 

Press  R/S  to  display  ENDR 

Press  R/S  again  to  display  f(X) 

Press  R/S  once  more  to  display  X 

To  repeat  program  press 
f CLREG  and  go  to  Step  5 


Key(s)  Output 


User  supplied 
function 
program 


2 . Sample  Problem 


2 

minimize  f(x)  = |x  - 16 1 
subject  to  x t Iq 

IQ  = 11,161,  RIOU  = 0.01,  NMAX  = 25. 

SOLUTION:  [3.99560,  4.00240]  in  17  function  evaluations 

Function  Value  = 0.00162 
Minimum  Point  x = 3.99980 


142 


User  Supplied  Proqram  Listin 


139 

21  16 

4 C 
A V 

I4C 

r:le 

36 

IS 

141 

Si 

53 

142 

1 

61 

143 

6 

86 

144 

- 

-4f 

145 

06$ 

16 

31 

146 

24 

147 

R S 

51 

Recalls  x from  Register  E. 


Function  value  is  left  in  the 
display  (x-reoister) . 


1.68  CS6n 
le.66  £$£>£ 
.61  6SE2 
25.66  6S62 
6SBE 
3.99568  »*« 
R<  $ 

4.60246  **» 
k 6 

6.8C162  »m 
R 5 

3.99988  »#* 


First  load  provided  proqram  (with  RUN  position) . 
Select  GTO  f e.  Move  switch  to  W/PRGM.  Enter  user 
supplied  proqram  for  function.  Move  switch  to  run 
position.  Enter  input  and  compute . 


b 


Minimize  f(x)  - 20  - x + ■■ 

(16-x) 


Subject  to  x ( IQ 

I0  * [10,15.9],  RIOU  - 0.01,  NMAX  = 25 

SOLUTION:  [14.99555,  15.00255]  in  15  function  evaluations 


Function  Value  = 6.00000 
Minimum  Point  x = 14.99823 


User  Supplied  Proqram  Lis tin 
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$LBie 

21  lb  15 

141? 

2 

32 

141 

3 

ec 

142 

RCLE 

36  15 

143 

- 

-45 

144 

i 

6! 

145 

b 

86 

14t 

RCLE 

' 3e  15 

147 

- 

-45 

14? 

l'X 

53 

149 

+ 

-55 

158 

RTH 

24 

1 51 

R'  S 

51 

Recall  x from  Recrister  E. 
Result  is  left  in  x-register 


18.83  iSBr! 
15.98  6666 
• Gl  £5Bl 
c5« 8^  fc$£. 
6 SEE 

14.995 55  **» 
fi/S 

15.P0255  »** 

R'  S 

6.88380  mi 

R'S 

14.59*23  *»* 


First  load  provided  program  with  switch  in  RUN  position 
Select  GTO  f e.  Move  switch  to  W/PRGM.  Enter  user 
supplied  program  for  function.  Move  switch  to  RUN 
position.  Enter  input  and  compute. 


minimize  f(x)  = j + sin(^)  (x  in  radians) 


subject  to  x € I 


I = 10,10],  RIOU  = 0.01,  NMAX  = 25, 


SOLUTION:  [2.79094,  2.79827]  in  16  function  evaluations 

Function  value  = 0.44890  j local  minimum 
Minimum  point  x = 2.79373  ( global) 


User  Supplied  Program  List in 


13? 

*LBL* 

21  16  15 

148 

RCLE 

36  15 

141 

62 

142 

T 

-24 

143 

Eh  T? 

-7  f 

a 

144 

EH7t 

-21 

145 

Pi 

16-24 

148 

X 

-35 

147 

SIN 

41 

148 

+ 

-55 

149 

FTH 

24 

156 

R'S 

51 

Recalls  R5  for  x utilization  and 
stores  function  value  in  R7 
(GSB  1) 


<3.80  fStn 
18.00  £SS£ 
.01  tSfC 
25.88  csb: 
£££>£ 
2.79694  *♦? 

R.i 

2.79S27  **, 

R.i 

8.44896  *»* 

*.-$ 

2.79273  *»* 


NOTE:  First  load  provided  program  with  switch  in  RUN  position 

Select  GTO  f e.  Move  switch  to  W/PRGM.  Enter  user 
supplied  program  for  function.  Move  switch  to  RUN 
position.  Enter  input  and  compute. 


3.  D*.ogram  Storage 

Allocation 

Registers : 

RO  : 

function  counter 

SO: 

RA: 

Rl: 

ENDL 

SI: 

RB: 

R2 : 

ENDR 

S2 : 

BC: 

R3 : 

RIOU 

S3: 

RD: 

R4  : 

NMAX 

S4: 

RE: 

Current  X 

R5  : 

X1 

S5: 

RI : 

4 (decrement) 

R6: 

X2 

S6 : 

R7 : 

F1 

S7 : 

R8 : 

F2 

S8: 

R9 : 

not  used 

S9  : 

NOTE:  User  supplied  function  can  utilize  sixteen  registers, 

Initial  Flag  Status  and  Use: 


0: 

OFF,  Unused 

2: 

OFF, 

Unused 

1: 

OFF,  Unused 

3: 

OFF, 

Unused 

User  Control  Keys 

A: 

Left  endpoint  (ENDL) 

a : 

B: 

Right  endpoint  (ENDR) 

b : 

C: 

Required  interval  of 
uncertainty  (RIOU) 

c : 

D: 

Maximum  function 
iterations  (NMAX) 

d: 

E: 

Compute 

e : 

User 

defined 

661 

*LBU 

21  11 

Input : 

002 

£701 

35  01 

Left  endpoint 

003 

RTH 

24 

Right  endpoint 

004 

*L  BLB 

21  12 

Required  interval  of  uncertainty 

005 

ST  02 

35  02 

Maximum  number  of  function  evaluations 

000 

RTH 

24 

00 7 

tLBLC 

21  13 

(ENDL,  ENDR,  RIOU,  NMAX) 

008 

ST03 

35  63 

009 

RTH 

24 

010 

tLBLD 

21  14 

011 

ST04 

35  64 

OH 

3 

63 

013 

014 

EHTt 

5 

-21 

65 

Calculate  Sigma 

015 

016' 

01? 

IX 

2 

54 

-45 

62 

a . <3 : ^ 

2 

BIS 

T 

-24 

019 

RCL2 

36  02 

020 

021 

RCL1 

36  01 
-45 

Calculate  initial 

022 

X 

-35 

823 

024 

RCL 1 
+ 

36  61 
-55 

X1  = ENDL  + a (ENDR  - ENDL) 

025 

ST05 

35  05 

028 

027 

0 

sroe 

66 
35  66 

Initialize  counters 

828 

028 

4 

STOI 

64 

35  43 

Set  up  display 

030 

RCL  4 

36  64 

Require  radian  calculations 

031 

DSF5 

-63  05 

032 

RHD 

16-22 

033' 

RTH 

24 

034 

tLBLE 

21  15 

035 

RCL  5 

36  65 

036 

RCL1 

36  61 

03? 

038 

RCL2 

-45 
36  62 

Calculate  initial  X^ 

039 

XZY 

-41 

040 

- . 

-45 

041 

ST06 

35  66 

, * 

042 

043 

044 

SSB1 
6SB8 
GSB  2 

23  61 
23  68 
23  62 

% 

Obtain  initial  function  values 

(F^  and  F^)  from  initial  X^  and  X^. 

045 

GSB  8 

23  62 

046 

»IPU 

21 

16  1 1 

047 

RCL8 

36  06 

846 

FCL7 

36  07 

84? 

X'  Y9 

16-34 

85  8 

6703 

22  0: 

051 

FCL6 

36  66 

053 

ST02 

35  02 

853 

FCL1 

36  01 

854 

- 

-45 

055 

RCL3 

36  03 

056 

x:y 

-41 

857 

X£Y? 

16-35 

856 

6T04 

22  04 

85? 

R CL  5 

36  05 

868 

ST06 

35  06 

061 

R CL  7 

36  or 

862 

ST  08 

35  06 

863 

CS  P6 

23  06 

064 

ST05 

35  05 

065 

RCL4 

36  04 

666 

RCLO 

06  00 

867 

X>V?' 

16-34 

868 

croc 

22  00 

86? 

0881 

23  0: 

878 

6S8S 

23  06 

■ 871 

6T0c. 

22 

16  11 

872 

*LBL3 

21  03 

873 

RCL2 

36  62 

874 

RCL5 

36  05 

875 

ST  01 

35  01 

876 

- 

-45 

877 

RCL3 

36  05 

878 

x:y 

-41 

87S 

X4Y9 

16-35 

086 

6T04 

22  04 

881 

RCL6 

36  06 

8S2 

STu5 

35  05 

883 

RCL8 

36  06 

084 

sto: 

35  or 

885 

6886 

23  06 

866 

ST  06 

35  06 

887 

RfL4 

36  04 

888 

RCLO 

36  00 

869 

X>  V° 

16-34 

896 

e roc 

22  60 

Determine  larqer  function  value.  if  F^ 
larger  go  to  branch  two,  otherwise  go  to 
branch  one. 


X.,  becomes  riqht  endpoint  of  interval  of 
uncertainty 


Compare  10b  RIOU? 

If  yes,  print  out  final  results 


Old  X^  becomes  X., 
Old  F^  becomes  F^ 


Determines  new  X **  (X^)  and  store. 


lias  NMAX  been  exceeded?  If  yes,  print  error. 


Calculate  Fj  arui  return  to  compare  new 
values  of  F'i  and  F2. 

Xj  becomes  left  endpoint  of  interval 
of  uncertainty 


Compare  IOU  Rlou? 

If  yes,  print  out  final  results 


Old  X,  becomes  X^ 
Old  F,  becomes  F^ 


Determines  new  X (X,)  and  store. 

lias  NMAX  been  exceeded? 

If  yes,  print  error. 


Calculate  F,  and  return  to  compare  new 
values  of  F.  and  F. 


Increments  function  counter  for 
of  NMAX  e xce e de d . 

determination 

Determines  new  X (X^  or  X, 

depending 

on  which  branch  subroutine  called  from) 

Determines  which  function  value 

with 

correspond inq  KNDL,  KNDR,  and 

X should 

be  printed.  (Smaller  function 
used. ) 

value  used.) 

Sets  up  stack  for  printout  if 

small  value. 

F1 

Sets  up  stack  for  printout  if 

smaller  value 

F2 

Prints  out  final  results. 

Fj  Rout ine 

F,  Routine 

Use  defined  function 


1 


E . Mathematical  Analysis 

For  each  X the  function  f(X)  is  evaluated.  By 
comparing  two  function  values  and,  F2  the  interval 

of  uncertainty  can  be  reduced  as  follows: 


[ 

i 


f 


6 


0 


jimimi'i  -■  ♦ ■ — jj—  x 


The  placement  of  the  X's  is  determined  by  the  Golden  Section 
Search  Technique  utilizing  the  "Golden  Ratio."  Placement  of 
the  first  X (X^)  determines  the  placement  of  all  other  X's 
since  all  remaining  points  are  placed  symmetrically  with  respect 
to  each  point  remaining  in  successive  intervals  of  uncertainty. 

In  order  to  ensure  that  each  IOU  length  is  predictively 
independent  of  the  function  f(x),  X^  and  X2  are  placed 
symmetrically  in  the  IOU.  When  the  IOU  is  reduced,  the  new 
shorter  IOU  still  contains  the  best  point  thus  far  achieved, 
so  selecting  a new  X point  will  allow  further  reduction. 

(New  X for  each  iteration.) 

Golden  Section  Search  selects  X^  to  satisfy: 


_ 


The  relative  position  of  the  X points  in  the  remaining 
IOU  is  the  same  at  each  iteration. 

Explanation : 

oL0 

I b 1 Iq  with  X. 

a Xx  b 1 

L0 

~ 1q  since  no  reduction  with  one  point.  Relative  position 
of  Xx  in  ^ is  glq/Lq  = ° . Now  place  X2  symmetrically 


OL, 


( 1- 2o ) L 


0 


oL, 


X0  X1 


Suppose  f(X2)  < f(X^)  and  reduce  IOU  accordingly 
(l~2o) Ln 

I 1 j2 

a X2  b 

v ' 

L2  = ( 1 -a ) 

Relative  position  of  X2  is 


15  3 


r 


(1  - 2o)L- 


1 - 2a 


(1  - a)  Lq  1 - o 


= ~ -2-l  ==>  ("relative  position  same  at  each  iteration") 

l-o 


3 + & 


a*  - 3o  + 1 = 0 ==>  o = 


use  root  o = 


3 - 


Choose  = ob  + (l-o)a  = a + o(b-a). 


Example:  n = 4, 

f ( x)  = 

(X  - . 303) 

2,  i0  = to,i] 

• 

K 

! 

W 

XK 

f(xK) 

lk 

0 

— 

— 

— 

[0,1] 

1 

1 

[0,1] 

. 382 

( .079)  2 

[0,1] 

1 

2 

[0,1] 

.618 

( . 31 5 ) 2 

[0,  .618] 

.618  = 

(l-o) 

3 

[0, .618] 

.236 

( .067)  2 

[0, .382] 

.382  = 

(l-o) 

4 

[0,  . 382] 

.146 

(.157)  2 

[.146, .382] 

.2  36  = 

(l-o) 
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